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INTRODUCTION 


A.  Historical  Background 

The  history  of  Fourier  transform  techniques  in  optic  dates  back  to 
the  experiments  of  Abbe  (1873)  [1],  [2]  and  later,  Porter  (1903)  [3].  These 
experiments  dealt  with  manipulating  the  spatial  frequency  content  of  a  scene. 
As  an  example  of  the  results  obtained  in  these  early  investigations,  consider 
an  input  scene  consisting  of  a  regular  array  of  squares  as  shown  in  Figure  1. 
If  this  is  illuminated  with  reasonably  monochromatic,  collimated  light  and 
focused,  the  spatial  frequency  spectrum  of  Figure  1  results  in  the  focal 
plane.  Abbe  and  Porter  were  the  first  to  discover  that  an  obstacle  placed  in 
this  plane  could  dramatically  alter  the  reimaged  scene.  A  small  slit  placed 
in  the  focal  plane  in  the  horizontal  direction  was  shown  to  cause  the  omission 
of  the  horizontal  components  of  the  resulting  reimaged  scene.  The  same  slit, 
allowing  only  the  vertical  spectra  to  pass,  produced  a  horizontal  bar  pattern, 
filtering  out  the  vertical  pattern.  This  is  illustrated  in  Figures  2  and  3. 
Abbe  and  Porter  also  observed  that  when  a  small  stop  is  placed  in  the  focal 
plane  blocking  the  lowest  order  (central)  spatial  frequency  component,  a 
contrast  reversal  took  place  in  the  resulting  image.  These  very  basic  experi¬ 
ments  had  far-reaching  effects  in  the  development  of  phase  contrast  microscopy 
as  proposed  by  Fritz  Zernike  in  1935  [4].  The  desire  was  to  construct  a 
microscope  so  that  the  intensity  reaching  the  viewer  is  linearly  related  to 
the  phase  change  undergone  by  the  light  in  passing  through  the  subject  under 
study.  This  was  to  prove  very  useful  in  the  study  of  nearly  transparent 
subjects.  Zernike  achieved  this  result  by  placing  a  “phase  changing  plate”  ir. 
the  focal  plane  of  the  microscope.  This  plate  can  be  a  small  dot  on  a 
transparent  medium  with  the  dot  centered  on  the  focal  point.  The  thickness  of 
the  dot  and  its  index  of  refraction  can  be  chosen  so  that  it  retards  the  phase 
of  the  focused  illumination  by  ir/2  (or  3ir/2)  relative  to  the  phase  retardation 
of  the  light  diffracted  from  the  subject.  Thus,  the  image  intensity  becomes 
related  linearly  to  the  phase  shift  caused  by  the  subject.  It  is  by  this 
method  that  we  are  able  to  view  a  phase  shift  as  a  change  in  intensity. 

In  1948,  Dennis  Gabor  published  a  paper  setting  down  the  basic  for¬ 
mulation  of  what  he  called  holography  [5].  He  reasoned  that  both  amplitude 
and  phase  information  could  be  recorded  if  the  illuminating  coherent  light 
could  be  interfered  with  the  illumination  carrying  the  image  information.  The 
basic  arrangement  for  making  a  Gabor  hologram  is  shown  in  Figure  4.  In  making 
a  Gabor  hologram,  the  object  must  consist  of  a  semitransparent  scene  on  a 
transparent  background,  a  limitation  that  was  not  removed  for  14  years.  The 
transmitted  light  then  consisted  of  a  scattered  portion  produced  by  the  image 
and  an  unscattered  part  transmitted  by  the  transparent  part  of  the 
transparency.  The  two  then  interfered  at  the  film.  This  mechanism  is 
discussed  in  Section  II.  To  view  the  reconstructed  image  one  had  only  to 
reilluminate  the  developed  hologram  transparency  as  shown  in  Figure  5.  The 
Gabor  hologram  did  possess  two  real  weaknesses.  As  mentioned  before  the 
object  had  to  be  quite  transparent.  This  severely  limited  the  types  of 
objects  that  could  be  used.  A  second  problem  which  is  shown  in  Figure  5  is 
the  generation  of  a  so-called  twin  image.  The  illuminated  hologram  not  only 
transmitted  the  Incident  light  but  also  reflected  a  sizeable  portion  which 
produced  an  image  on  the  other  side  of  the  hologram  and  the  two  could  not  be 
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separated.  If  one  focused  on  the  real  image,  the  out-of-focus  virtual  image 
was  present  and  vice  versa.  For  all  of  its  limitations  the  Gabor  hologram  did 
set  the  foundation  for  the  development  of  modern  holography. 

It  was  also  at  about  this  time  that  the  concept  of  a  matched  filter 
was  introduced  by  North  in  1943  [6] ,  and  applied  to  signal  detection  by 
Peterson,  Birdsall,  and  Fox  in  1954  [7].  A  matched  filter  is  described  by 
these  authors  to  be  the  optimum  linear  filter  for  signal  detection  in  white 
noise.  Of  course,  these  early  matched  filters  were  used  for  the  detection  of 
electrical  signals  embedded  in  noise  but  the  concepts  carried  over  into  opti¬ 
cal  processing. 

Perhaps  the  most  dramatic  successes  in  optical  spatial  filtering 
came  in  the  area  of  photograph  improvements.  Marechal  at  the  Institut 
d'Optique,  Universite  de  Paris  in  the  early  1950's  is  probably  the  most 
noteworthy  [8],  He  believed  that  defects  in  photographs  arose  from  problems 
in  the  ability  of  the  photographic  system  to  faithfully  transfer  the  spatial 
frequency  content  of  the  original  scene.  He  attempted  to  insert  appropriate 
attenuating  and  phase  shifting  plates  in  the  focal  plane  in  order  to  develop 
filters  which  would  remove  some  of  the  defects.  Probably  Marechal's  most 
important  discovery  was  that  small  details  in  the  image  could  be  enchanced  if 
the  low  spatial  frequency  components  could  be  blocked  out  or  attenuated.  The 
reverse  of  this  was  also  observed  in  the  filtering  of  the  periodic  structure 
associated  with  the  halftone  printing  process.  This  produced  a  continuous 
image  rather  than  an  image  obviously  made  up  of  black  dots  on  a  white 
background.  Marechal  also  made  a  considerable  contribution  to  the  removal  of 
image  blur  by  using  an  absorbing  plate  and  a  dephasing  plate  in  the  focal 
plane  of  the  filtering  system. 

In  the  decade  of  the  1950's  the  similarities  between  the  fields  of 
optics  and  electrical  engineering  (particularly  in  communications)  became 
obvious  and  a  great  deal  of  exchange  began  to  take  place.  The  problems  asso¬ 
ciated  with  filtering  and  detection  shared  common  approaches  to  solution.  A 
paper  published  in  the  Journal  of  the  Optical  Society  of  America  by  a  com¬ 
munications  theorist,  Peter  Elias  and  others  in  1953  [9]  prompted  the  initial 
exchange  of  ideas.  The  combining  of  the  fields  was  perhaps  best  provided  by  a 
physicist,  E.  L.  O'Neill  in  1956  [10]  and  by  a  symposium  in  1960  entitled 
"Communication  and  Information  Theory  Aspects  of  Modern  Optics"  [11].  At  pre¬ 
sent  the  fields  are  quite  closely  intertwined  and  one  will  find  articles  in 
electrical  and  communication  engineering  journals  written  by  opticists  and 
articles  in  optical  journals  written  by  engineers.  It  all  comes  down  to  the 
methods  by  which  information  is  transferred  from  one  point  to  another  - 
whether  it  be  by  electrical  or  optical  means. 

In  1960,  the  advances  ir.  optics  and  communication  engineering  began 
to  accelerate  rapidly.  The  laser  was  developed  [12]  which  gave  the  fields  a 
new  source  of  high  intensity  coherent  light  to  use.  Radar  signal  processing 
using  coherent  filtering  got  underway  at  the  University  of  Michigan  Radar 
Laboratory.  The  publication  of  a  paper  entitled,  "Optical  Data  Processing  and 
Filtering  Systems"  in  1960  by  this  group  [13]  under  the  direction  of  L.  J. 
Cutrona,  was  to  be  the  first  of  many  to  draw  the  attention  of  both  engineers 


and  physicists.  The  military's  interest  in  the  area  of  optical  data  pro¬ 
cessing  began  at  about  this  point  in  time.  The  Michigan  group  was  funded  to 
investigate  the  processing  of  data  collected  by  “synthetic  array"  antennas 
aboard  reconnaissance  planes  [14].  In  order  to  use  microwave  radar  and  get 
reasonable  resolution,  the  size  of  the  antenna  would  be  much  too  large  to  be 
flown  aboard  a  surveillance  plane.  The  resolution  at  a  range  R  goes  approxi¬ 
mately  like  XR./D,  where  X  is  the  wavelenght  and  D  is  the  linear  size  of  the 
antenna.  For  a  resolution  of  one  meter,  a  plane  flying  at  an  altitude  of  3  km 
would  need  an  antenna  450  m  long  if  the  radiation  used  had  a  wavelength  of  15 
cm.  In  order  to  solve  this  problem,  a  system  was  developed  which  used  pulses 
of  radiation  so  that  the  information  obtained  is  a  periodic  sampling  of  the 
terrain.  The  task  of  the  Michigan  group  wau  to  turn  this  sampled  information 
into  a  real  image  which  they  did  using  the  data  obtained  to  form  a  moving  film 
transparecny ,  then  demodulating  it  using  an  optical  system  of  their  own  device. 

This  brief  survey  is  by  no  means  complete.  Optical  processing  has 
been  used  successfully  in  many  other  fields  such  as  Fourier  spectroscopy  and 
seismic  wave  analysis,  both  of  which  have  benefited  greatly  from  the  basic 
analysis  methods  discussed  here. 

B.  Recent  Investigations 

This  section  must  surely  be  in  with  the  application  of  Gabor's 
theories  and  experiments  in  holography  by  Leith  and  Upatnieks  of  the 
University  of  Michigan's  Radar  Laboratory  in  1962  [15].  Although  this  is 
after  the  development  of  the  laser  in  1960,  the  first  investigations  were  done 
without  a  laser.  The  Leith-Upatnieks  hologram  was  a  vast  improvement  on 
Gabor's  in  that  it  removed  two  of  the  most  serious  limitations  of  the  Gabor 
hologram:  the  requirement  of  a  transparent  object  and  the  production  of  an 
inseparable  image  twin.  The  basic  experimental  arrangement  for  making  a 
Leith-Upatnieks  hologram  is  given  in  Figure  6.  The  noticeable  difference  be¬ 
tween  the  Gabor  hologram  and  the  Leith-Uptanieks  hologram  is  the  separate 
reference  beam.  In  Figure  6  a  semitransparent  object  is  used,  but  even  that 
is  not  necessary.  An  opaque  object  or  even  a  real  solid  model  may  be  used  and 
the  light  collected  by  reflection  instead  of  transmission.  Light  from  the 
object  intefers  with  the  reference  beam  at  the  location  of  the  film  plate. 

After  developing  the  film  the  recreated  image  may  be  viewed  by  reillumination 
with  the  source.  The  image  occurs  in  focus  off  axis  along  the  path  of  the 
original  reference  beam.  The  mathematical  details  of  the  analysis  of  this 
type  system  are  left  to  Section  II.  It  is  also  interesting  to  note  that  the 
first  hologram  using  a  laser  as  the  source  was  produced  by  this  group  [16]. 

One  year  later,  in  1963,  A.  B.  VanderLugt  of  the  University  of 
Michigan's  Radar  Laboratory  proposed  and  demonstrated  a  new  technique  for 
synthesizing  frequency  plane  masks  for  coherent  processors  [17].  The  filters 
constructed  using  this  technique,  like  holograms,  have  the  unique  property  of 
being  able  to  control  both  the  amplitude  and  phase  of  a  transfer  function. 

The  technique  used  in  making  these  filters  is  almost  identical  to  that  used  in 
making  an  ordinary  hologram  except  that  the  input  scene  is  Fourier-transformed 
by  a  lens  before  it  is  allowed  to  interfere  with  a  reference  beam  at  a 
holographic  film  plate.  Properly  stated,  filters  made  using  this  technique 


7 


should  be  called  holographic  Fourier  transform  matched  filters.  A  sketch  of 
the  experimental  arrangement  used  in  making  this  type  of  matched  filter  is 
given  in  Figure  7. 

After  the  exposed  film  is  developed,  it  is  placed  in  an  optical 
system  as  sketched  in  Figure  8.  If  the  input  scene  is  the  same  one  used  to 
record  the  matched  filter,  a  correlation  intensity  will  be  recorded  by  the 
detector  along  the  axis  of  the  original  reference  beam.  The  mathematical 
details  of  this  phenomena  are  given  in  Section  II.  The  important  point  to  be 
made  here  is  that  using  the  previously  described  experimental  arrangement,  an 
optical  identification  device  was  first  demonstrated.  The  VanderLugt  tech¬ 
nique  for  making  matched  filters  is  the  basis  for  virtually  all  the  work  done 
in  this  area  today. 

In  the  time  from  1965  to  1970,  many  authors  investigated  a  variety 
of  topics  dealing  with  optical  matched  filtering.  One  of  the  most  pressing 
problems  was  that  of  multiple  recognition.  The  VanderLugt  arrangement  just 
described  was  limited  to  the  recognition  of  one  input  scene  per  filter. 
LaMachia  and  White  proposed  a  method  of  multiplexing  using  angularly  separated 
reference  beams  [18].  VanderLugt  suggested  a  similar  method  of  superimposing 
a  number  of  filters  with  each  filter  on  a  different  carrier  [19].  Leith 
devised  a  technique  of  coding  the  superimposed  filters  by  rotating  the 
holographic  plate  for  each  exposure  [20],  Groh  used  point  holograms  to  pro¬ 
duce  spatially  separated  filters  with  different  carriers  [21].  Later,  Shi 
proposed  color  multiplexing  using  a  multiple  wavelength  source  [22].  In  1972, 
Grumet  developed  an  NxN  array  of  holographic  lenses  which  were  used  to  Fourier 
transform  a  scene  in  a  VanderLugt  type  optical  correlator  [23]. 

From  1972  to  1977,  considerable  effort  was  placed  on  determining  the 
usefulness  of  optical  matched  filers  in  object  recognition  and  detection.  In 
1972,  the  Harris  Electro-Optics  Center  of  Radiation  produced,  under  contract 
to  the  US  Army  Missile  Command,  a  report  discussing  the  possible  use  of  opti¬ 
cal  matched  filters  for  terminal  missile  guidance  [24],  The  technical  aspect 
of  this  research  was  under  the  direction  of  A.  B.  VanderLugt.  The  investiga¬ 
tion  concerned  incoherent  (mask  matching)  correlation  and  coherent  (holo¬ 
graphic  Fourier  transform)  correlation.  Conclusions  drawn  from  the  research 
include:  Incoherent  system  suffers  two  major  drawbacks;  (1)  diffraction 

effects  and  the  presence  of  a  large,  often  nonuniform  bias  in  the  output,  and 
(2)  the  coeherent  system  was  shown  to  have  little  of  either  problem.  Experi¬ 
ments  were  performed  which  showed  that  the  correlation  intensity  increases 
with  the  resolution,  spatial  frequency,  of  the  input  scene  for  a  coherent  or 
incoherent  system.  Rotation  and  scale  change  experiments  were  also  done  which 
showed  that  the  coherent  correlation  intensity  fell  by  13  dB  when  the  input 
scene  was  rotated  by  about  one  degree  and  the  intensity  fell  by  20  dB  with  a 
1  percent  change  in  scale.  The  input  scene  used  was  an  urban  scene  containing 
moderately  high  spatial  frequencies.  Results  of  several  coherent  and  inco¬ 
herent  matched  filter  experiments  were  given  which  show  the  definitive  advan¬ 
tage  of  the  coherent  system,  especially  for  high  spatial  frequency  input 
scenes.  This  report  was  also  one  of  the  first  to  mention  the  need  for  a  real 
time  correlator  without  specially  prepared  transparencies  for  input  scenes.  A 
photo-plastic  device  developed  by  Harris  Electro-Optics  was  tested  and  found 
to  be  moderately  useful  but  slow  in  reacting  (about  5  seconds). 
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In  1977,  Casasent  and  Furman  published  one  of  their  many  papers 
dealing  with  optical  correlation  [25].  This  particular  paper  is  of  interest 
because  they  investigated  the  effects  of  rotation,  scale  and  filter  shifts  on 
the  correlation  and  obtained  results  similar  to  those  of  Vander  Lugt's  group, 
just  discussed.  The  low  resolution  scene  used  was  an  aerial  map  of 
Huntsville,  AL  having  a  resolution  of  details  on  the  map  of  ~  30  pm.  The 
magnification  of  the  map  Itself  was  not  given. 

This  map  contained  both  urban  (high  spatial  frequency)  and  rural 
(low  spatial  frequency)  scenery.  The  input  scenes  were  35  mm  slides  (Kodak 
high  speed  holographic  50-153  film).  The  dimensions  of  the  photographs  were 
not  given  so  an  estimate  of  the  spatial  frequency  content  of  the  slide  itself 
is  not  possible,  although  it  appears  to  be  near  that  of  the  Harris 
investigation. 

In  1978,  Grumman  Aerospace  concluded  an  investigation  into  using 
optical  Fourier  transform  matched  filters  for  aerial  reconnaissance  film 
screening.  Kenneth  Leib  was  the  principal  investigator.  This  team  developed 
the  first  application  of  a  multifocus  holographic  element  to  matched 
filtering.  In  the  1978  article  discussing  this  work,  the  investigators 
reported  on  a  holographic  element  capable  of  producing  a  2X2  array  of  Fourier 
transforms  addressable  with  35  mm  slides  (Kodak  SO-173  film)  [26].  The  pri¬ 
mary  input  scene  of  interest  was  a  scale  model  of  a  US  Army  M-60  tank,  but 
later  some  terrain  imagery  was  used.  The  noteworthy  facet  of  this  investiga¬ 
tion  was  the  attempt  at  multiple  exposures  of  the  input  scene  while  using  the 
2X2  array  (Fourier  transform  image  overlaps).  They  found  that  no  more  than 
four  Fourier  transforms  could  be  stored  at  the  same  spot  without  significant 
loss,  84  percent,  in  correlation  when  compared  to  the  correlation  from  a 
single  stored  scene.  This  study  was  done  using  the  terrain  scene  of  moderate 
spatial  frequency  content.  A  study  of  the  rotational  dependence  of  the  corre¬ 
lation  was  done  which  showed  that  for  the  low  spatial  frequency  input  (M-60 
tank) ,  the  scene  could  be  rotated  by  more  than  7  degrees  before  the  correla¬ 
tion  dropped  by  50  percent.  This  is  in  contrast  to  the  work  by  Harris  and 
Casasent  where  only  tenths  of  degrees  were  needed  for  the  same  decrease.  The 
wide  difference  in  spatial  frequencies  of  the  input  scenes  may  explain  the 
difference  in  results. 

In  1978,  another  important  modification  to  the  optical  correlator 
was  made  by  the  addition  of  an  lncoherent-to-coherent  image  converter,  liquid 
crystal  light  valve  (LCLV) ,  first  described  by  Hughes  in  1975.  This  device 
allowed  the  use  of  real  input  scenes  such  as  those  viewed  by  the  eye.  The 
making  of  a  transparency  was  no  longer  necessary.  The  device  was  reasonably 
fast-reacting,  allowing  25  ms  "real  time"  correlations  for  the  first  time  [28]. 
This  device  is  described  in  Section  III. 

This  experiment  also  demonstrated  for  the  first  time,  millisecond 
tracking  of  a  target.  An  interesting  facet  of  optical  correlation  is  its 
translational  invariance,  that  is,  the  input  scene  can  be  moved  without  losing 
the  correlation;  in  fact  the  correlation  spot  will  move  along  with  the  input 
scene.  This  was  demonstrated  quite  vividly  in  a  paper  by  Guenther  [29]  and  in 
a  video  tape  produced  by  Farr  [30].  The  real  time  correlator  could 
distinguish  between  two  similar  vehicles  and  track  a  vehicle  as  it  moved 
through  the  field  of  view. 


Christensen,  in  1979,  showed  that  up  to  eight  matched  filters  could 
be  stored  at  a  single  focal  point  using  multiple  exposure  techniques  [31]. 

This  was  done  with  the  aid  of  an  LCLV  making  it  real  time  recognition.  Leib, 
of  Grumman  Aerospace,  was  able  to  store  4  filters  using  a  holographic  lens 
array  [32].  This  was  not  done  using  an  LCLV,  but  did  demonstrate  the 
parallel  -  addressing  ability  of  such  a  system  -  more  than  one  matched  filter 
could  be  addressed  at  a  time. 

More  recently,  Liu  and  Duthie  have  proposed  a  multiplexing  method 
incorporating  a  special  contact  screen  in  a  real  time  matched  filter  correla¬ 
tor  [33].  The  screen  used  in  that  investigation  is  described  in  Section  II. 
Essentially  the  screen  split  an  input  scene  into  an  NxN  array  of  identical 
scenes,  which  were  Fourier-transformed  using  a  conventional  lens  in  a  Vander 
Lugt  type  correlator.  Further  developments  in  this  approach  finally  led  to  a 
single  holographic  element  which  could  both  split  the  scene  into  an  array  and 
Fourier-transform  it.  Results  and  a  complete  description  are  given  in  a  paper 
by  Gregory  and  Liu  [34]. 

One  of  the  newest  fields  is  that  of  white  light  optical  processing. 

A  promising  techique  by  Yu  employs  a  high  diffraction  efficiency  grating  in 
contact  with  the  input  image  transparency  along  with  an  achromatic  transform 
lens  [35].  The  diffraction  grating  separates  the  white  light  into  its  color 
components  and  the  lens  transforms  each  component. 

The  problem  of  recognition  has  been  addressed  by  computer  experts  in 
some  detail.  At  present  digital  processes  of  recognition  are  being  used  for 
such  diverse  applications  as  lung  tumor  detection  and  chromosome  karyotyping 
[36].  This  technique  requires  some  highly  sophisticated  computer  hardware  and 
software  and  must  be  tailored  for  the  specific  recognition  problem.  In  addi¬ 
tion  to  being  quite  slow,  these  systems  have  a  fairly  high  error  associated 
with  them.  It  is  hoped  that  many  of  the  recognition  problems  now  handled 
digitally  can  be  done  optically  if  enough  matched  filter  memory  can  be  stored 
and  accessed. 

The  military  has  contributed  greatly  in  funding  for  optical  recogni¬ 
tion  research.  The  desire  is  to  be  able  to  distinguish  between  friend  or  foe 
without  having  to  place  personnel  in  particularly  dangerous  environments.  A 
device  which  could  recognize  a  military  vehicle,  ship,  or  plane  and  recognize 
it  immediately  would  have  particular  application  in  unmanned  reconnaissance 
missions.  In  terminal  missile  guidance,  a  system  such  as  this  would  allow  the 
programming  of  particular  highly  localized  targets  into  the  memory  of  the 
guidance  system,  thereby  making  pinpoint  attacks  possible.  The  tracking  capa¬ 
bility  of  such  a  system  would  allow  real  time  identification  to  be  made  on 
moving  targets  as  well. 

The  goal  of  the  present  investigation  is  to  bring  together  many  of 
the  topics  just  discussed,  in  combination  with  a  new  holographic  element  pro¬ 
duced  by  Liu  [33],  in  an  attempt  to  investigate  a  large  memory  optical 
recognition  system  operating  in  real  time.  The  memory  consists  of  an  array  of 
holographic  Fourier  transform  matched  filters,  addressable  in  series  or 
parallel. 
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The  basic  holographic  equations  are  given  in  Section  II  along  with 
simple  models  of  the  correlation  of  scenes  that  have  been  tilted,  rotated,  or 
magnified.  An  outline  of  the  basic  principles  of  the  Liu  holographic  element 
is  also  given. 

The  experimental  procedures  followed  in  producing  the  matched  filter 
arrays  are  given  in  Section  III.  The  liquid  crystal  light  valve  (LCLV)  used 
in  this  investigation  is  described  and  its  operating  parameters  given  along 
with  a  list  of  the  other  basic  equipment  used. 

The  experimental  results  given  in  Section  IV  range  from  basic  tests 
of  the  holographic  element  and  the  LCLV,  to  the  storing  and  accessing  of  more 
than  seventy  different  matched  filters.  The  present  state  of  the  system 
allows  the  parallel-addressing  of  25  different  filters  or  the  series-addressing 
of  up  to  a  demonstrated  73  filters,  in  real  time,  with  good  correlations 
resulting. 


II.  THEORY 


A.  The  Fourier  Transform 

It  is  well  known  that  any  function  f(F)  that  represents  a  real  physi 
cal  phenomena  can  be  expanded  in  terms  of  orthonormal  functions  Un(F)  f37]: 


f(p)  =  2]  anUn(r)  •  O) 

n=  1 


where  the  coefficients  an  are  given  by: 
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n 
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U*  (O  f(F.)  dF  . 
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The  integration  limits  span  the  interval  of  the  function  f ( F) . 

As  an  example,  consider  the  orthogonal  functions  sine  and  cosine. 
The  expansion  of  a  function  in  terms  of  these  functions  is  called  a  Fourier 
series.  The  series  equivalent  of  Equation  (1)  is: 
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where  the  interval  in  F.  is  (-  a/2, a/2)  and. 
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m  =  0,  +  1,  +  2  on  the  interval  (-a/2,  a/2)  as  before.  An  expansion 

would  then  be  written: 


f(  r) 


=  JL  ^  e  *-  (  2  ,™  r  /  a 

m=-°° 


whe  re 


e-l(  2  ttj7ip  ’/a)  f  (?’  ) 


(7) 


(8) 


If  the  function  Is  well-behaved  and  at  least  piecewise  continuous,  let  the 
Interval  become  infinite  (a  -*■  m)  and  transform: 
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The  result  Is  the  well-known  Fourier  Integral: 
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Equations  (12)  and  (13)  are  the  mathematical  basis  for  a  wide 
variety  of  theories  related  to  signal  propagation,  electrical  or  optical. 

B.  Fourier  Transforming  Properties  of  Lenses 

The  analysis  of  how  a  simple  focusing  lens  performs  the  Fourier 
transform  of  a  scene  is  the  foundation  of  the  description  of  Fourier  transform 
matched  filtering.  Consider  the  sketch  of  a  spherical  lens  in  Figure  9.  The 
phase  delay  6,  experienced  by  an  incoming  spatially  and  temporally  coherent 
wave,  at  any  (x,y)  coordinate,  after  passing  through  the  lens,  is  calculated 
by  beginning  with  the  expression  for  the  velocity  of  light  Vg  in  a  medium 
having  an  index  of  refraction,  n:  [38]  (p.  77) 
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Substituting: 
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multiplying  both  sides  by  time,  t 
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thus : 
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Solving  for  <t>  : 


2nnA 


nAk 


(16) 


(17) 


(18) 


The  wave  also  travels  the  distance  Aq  -  A  in  open  air  (n  =  1).  Thus  the 
total  phase  delay  is:  [38]  (p.  77) 

<t>p  =  nAk  +  k(  Aq-A) 

=  kA(n-l)  +  kAg  .  (19) 

Now,  define  a  transform  of  the  form: 

i*T 

tj?(*,y)  -  e 

Substituting  ^  from  Equation  (19)  gives: 
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Figure  9.  Geometry  for  calculating  phase  delay  through  a  lens 


Thus  the  field  is  related  to  the  field  Uj>  by: 


U$(x,y)  =  tjj(x,y)  Ujj(x,y)  .  (21) 

The  problem  now  is  to  find  an  expression  for  the  thickness  function 
A(x,y).  Using  a  little  geometry  one  can  arrive  at  the  following  expression 
for  A(x,y)  [38]  (p.  77): 
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where  and  R2  are  the  radii  of  curvature  of  the  lens.  A  binomial  expansion 
(paraxial  approximation)  of  the  square  root  terms  gives: 
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and  similarly  for  the  term  involving  R2.  Substituting  this  into 
Equation  (22)  yields: 


tfx.y)  %  -  ^  ( I 


If  R2  =  -  =  R  (bi-convex  lens),  then: 

A(x,y)  *  Aq  -  x  +y. 


Note  that  for  any  thin  symmetric  lens: 


(24) 


(25) 


~  =  2(n-l)  |  (26) 

where  f  is  the  focal  length  of  the  lens.  Therefore,  the  final  form  of  the 
lens  transform  function  t  may  be  written  as: 


tjKx.y) 


a  ik  ,  2  2x 

iknAj,  -  2j  (x  +y  ) 


(27) 


To  apply  this  lens  transformation,  consider  the  typical  optical  system 
of  Figure  10.  The  input  plane  (F.,n)  contains  a  pattern  (scene)  stored  on  a 
transparency.  This  transparency  is  back-illuminated  with  coherent,  collimated 
light.  We  desire  to  know  what  the  image  looks  like  a  distance  f  away  from  the 
thin  spherical  lens. 


Between  the  (?.,rt)  plane  and  the  (x,y)  plane,  apply  Huygens  principle  [39]  (p. 
153): 


g(x,y)  =  C  JJ f(P,n)  elkrd?,dn 


where  g(x,y)  represents  the  field  at  the  entrance  to  the  lens,  s  denotes  the 
surface  integral,  and  C  is  an  arbitrary  complex  constant  related  to  the  ampli¬ 
tude  of  the  incident  field.  The  distance  r  may  be  written  as: 

r  -  [d2  +  (x-?)2  +  (y-n)2]^2  .  (29) 

Using  the  binomial  expansion  (paraxial  approximation)  gives: 

r  *  d  +  ~ -  [(x-O2  +  (y-n)2]  .  (30) 

Equation  (28)  may  be  rewritten  as: 


g(x,y)  =  C  f(?,,n)  *  h^F.n) 


where  C  is  again  a  complex  constant  and 


hl(F.,  n)  =  exp  £  f^<p.2+U2)] 


and  the  *  denotes  convolution.  In  general,  convolution  is  defined  as: 


f(z)  *  g(z) 


f(z  )  g( z-z ' )dz ' 


Upon  going  through  the  lens,  the  transform  t %  derived  earlier  may  be  applied 
to  give  the  field  as  the  wave  exits  the  lens: 


U*  -  g(x,y)  tj,(x,y) 


Traveling  from  the  lens  a  distance  f  to  the  screen  at  (xf,yj),  we  may  again 
apply  Huygens  principle.  Thus,  the  field  at  the  screen  can  be  written  as: 
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U(x,f,yf)  =  C  |  ^(P.n)  *  hj  (  r ,  n)J  t;(x,y)  j  *  h2(x,y) 


(35) 


where 

h2  =  exp  |^|y(x2+y2)J  (36) 

and  C' '  is  a  complex  constant.  Equation  (35)  may  be  rewritten  as: 


U(xf  ,yf) 


f(^,n) 


d?dn 


sj  s2 


(37) 


Sj  and  s2  denote  integration  over  the  (F.n)  plane  and  the  (x,y)  plane,  respec¬ 
tively.  The  term  y  is  given  by  [39]  (p.  155): 


Y  =  j  j^x-F.)2  +  (y-n)2J  +-i-[(xf-x)2  +  (yf-y)2  -  (x2+y2)J  .  (38) 

Using  the  binomial  expansion  again  produces: 

Y^sjfvE2  +  vx2  +  x2  -  2v*x  -  2xxf  +  vn2  +  vy2  +  y2 


where 


v 


(39) 


(40) 


After  a  great  deal  more  manipulation  [39]  (p.  156)  (not  involving  any  more 
approximations)  the  following  expression  for  the  field  at  the  screen  results: 


U(xf,yf) 


Ae*p  [^t  (^)  H+yf)] 
*  jU*f(F,ri)  exp[“FS'  (xf^+yf")] 


d?dn 


(41) 


Now  letting  the  distance  d  equal  f,  that  is,  placing  the  transparency  a  focal 
length  in  front  of  the  lens  (d  =  f  =>  v  =  1): 
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(42) 


U(xf.yf)  =  A  JJ  f(F,n)  exp  £  -p  (xfF.+yfn)J  dF.dn 


Comparing  Equation  (42)  with  Equation  (12),  we  see  that  the  true  Fourier 
transform  does  result.  Very  often  this  equation  will  be  written  in  the  form: 


U(xf,yf)  =  A 


// 


f(F,ri)  exp  |  -2iti  (Q^F.+OyO 


dF.dn 


(43) 


where  0X  and  Py  are  the  spatial  frequencies  of  the  input  scene  defined  by: 
[38]  (p.  61) 


xf 

Af  » 


(44) 


and  A  =  1 / i Xf  for  unit  amplitude  plane  wave  illumination  of  the  input  f(F,n). 
It  can  be  noticed  that  there  is  a  direct  correspondence  between  the  spatial 
measurements  and  frequencies  given  here  and  the  time  and  frequencies  used  in 
the  Fourier  analysis  of  time-varying  signals.  The  Fourier  transform  of  a 
scene  is  produced  by  a  lens  when  the  transparency  is  placed  a  focal  length  in 
front  of  the  lens  and  the  viewing  screen  a  focal  length  behind  the  screen. 
This  remarkable  result  forms  the  basis  for  all  Fourier  transform  matched 
filter  work  done  today. 


C.  Holography  and  Optical  Correlation 

Suppose  that,  in  the  previous  discussion,  the  viewing  screen  is 
replaced  by  a  photographic  film  plate  (such  as  Kodak  649F)  and  a  collimated 
reference  beam  (derived  from  the  original  source)  is  allowed  to  interfere  with 
the  Fourier  transformed  image,  thereby,  exposing  the  plate.  Figure  II  is  a 
schematic  representation  of  this.  The  reference  beam  may  be  represented  as  a 
plane  wave  having  amplitude  B  of  the  form:  [38]  (p.  172) 


-12  Tioryf 

Ur  =  B  e  (45) 

where  a  =  (sin  0)/A,  the  spatial  frequency  of  the  reference  beam.  The  total 
field  at  the  film  plate  is  then: 

UT  =  U(xf’yf>  +  Ur(yP  *  (46) 

Substituting  U(xf,yf)  from  Equation  (43)  and  Ur(yf)  from  Equation  (45)  yields: 

-i2xayf 

UT<xf *yf)  "  Yf  F  f  f^.Tl)}  +  B  e  (47) 

where  F  f  I  denotes  the  Fourier  transform.  The  intensity  of  these  combined 
beams  hitting  the  film  is: 

la  |  uT  |  2  -  ufuT  (48) 
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Figure  11.  Fourier  transform  hologra 


where  the  *  (super  star)  denotes  the  complex  conjugate.  Performing  this 
operation  we  arrive  at:  [381  (p.  172) 


I  =  const. 


o  7  R  i2  nay 

{  f(C.r)}  I'  +  +  —t  e  F{  f  ( r  . r  )  1 


-i2Tiay 

F*  {  f(r 


,n)  }  .  (4 


9) 


Suppose  for  a  moment  that  the  Fourier  transform  F  {f(^,ri)}  can  be  written  as 
an  amplitude  and  phase  product  such  as:  [38]  (p.  209) 


F  {  fU.n)  } 


A(F.,r.) 


n) 

e 


(50) 


much  the  same  way  as  the  Fourier  transform  of  a  time-varying  signal  is  often 
written  as  an  amplitude  and  phase  product.  After  substitution  and  a  little 
manipulation,  the  intensity  can  be  rewritten  as: 

r  1  2  2  2AB  ,  .  „  , 

I  =  const.  --  -  --  A  +  B  +  cos  (^  “  2iroty  ) 

L  \zr 

Notice  that  the  intensity  is  a  function  of  the  amplitude  and  phase. 
If  the  photographic  plate  is  exposed,  and  assuming  a  linear  response  by  the 
film,  then  the  developed  film  will  have  an  amplitude  transmission  proportional 
to  the  intensity  of  the  interfering  beams 


(51) 


t(xf,yf)  =  D  I(xf,yf) 


(52) 


where  t  is  the  transmission  and  D  is  the  film  response  constant.  This  is  the 
mechanism  by  which  both  amplitude  and  phase  information  may  be  stored  on  an 
intensity-sensitive  photographic  plate  [38]  (p.  210).  The  developed  film  is 
called  a  VanderLugt  or  holographic  Fourier  transform  matched  filter,  matched 
to  the  input  scene  (f(£,r))}.  If  this  developed  holographic  film  is  placed  in 
the  arrangement  sketched  in  Figure  12,  an  interesting  result  occurs.  The 
function  g(  t))  now  represents  a  test  input  scene,  maybe  different  from 
f(£.h).  The  field  F{g(^,p)}  falls  on  the  film  (again  assuming  coherent  plane 
wave,  unit  amplitude  illumination  of  g(^,t))  and  is  partially  transmitted. 

The  transmission  is  given  by  the  product  of  the  field  and  the  transmission  of 
the  plate  given  by  Equations  (52)  and  (49): 

U(trans)  =  F{g(x,y)}  •  | T { f 

where  the  constants  have  been  omitted. 

Equation  (49)): 


('.M|  |  2  (53) 

This  can  be  rewritten  as  (using 


■m 
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Figure  12.  Holographic  Fourier  transform  matched  filterin; 


where 


U(trans)  =  ~  +-^-T  F  r  G  +  —  F  G  e 
At  XJf J  Ki 


+  B  *_  -i2nayf 

xV  F  G  e 

(54) 

F|g(r,Tl)}  , 

(55) 

F  1  f  (  F, ,  n)  1  . 

(56) 

The  final  lens  L2  Fourier  transforms  U(trans)  so  that  the  field  In 
the  plane  (xc,yc)  Is  given  by:  [38]  (p.  175) 


Uc  =  B2g(xc,yc)  +  ["f(xc,yc)  *  f*  (-xc,-yc) 

B  *  f  L 

+—  f(xc,yc)  *  g(xc,yc)  *  6(xc,yc  +  aXf) 

+  ^  |f*(-xc  ,-yc)  *  g(xc,yc)  *  Mxc,yc  -nXf)J 
where  we  have  used  the  relation  [40] 

F  1-F  If  l*F  |G  1 


(57) 


(58) 


and  *  denotes  the  convolution.  The  6  function  comes  from  the  integration  of 
exp  (-12itayf).  Autocorrelation  is  defined  mathematically  to  be:  [39]  (p.  10) 


/ 


h(x) 


h*(x-a) 


dx 


—  OO 


(59 


The  final  term  of  Equation  (57)  may  be  written  as: 


ff 


g(  p  ,") 


f*  ()r-xc,n-yc  +  a\f)  dFdri 


(60) 


which  is  the  mathematical  definition  of  cross-correlation.  Note  that  the 
value  of  the  integral  is  zero  unless  yc  =  a\f.  Physically  this  means  that  the 
cross-correlation  of  g(F,,n)  and  f(F,n)  may  be  found  in  the  (xc,yc)  plane  but 
only  at  the  location  when  yc  =  a\f;  that  is,  when: 


yc 


sin9 

X 


(Xf)  , 


yc  -  f  sin  9 


(61) 

(62) 


The  cross-correlation  is  found  along  the  trajectory  of  the  original  reference 
beam  in  the  (xc,yc)  plane.  If  g(F,n)  =  f(P.,ri)  then  the  result  from  Equation 
(60)  is  called  the  auto-correlation. 


In  application,  the  correlation  (auto-  or  cross-)  is  a  measure  of  the 
similarity  of  two  scenes.  If  we  have  two  transparencies  of  a  scene,  and  lay 
one  over  the  other,  the  amount  of  overlap  of  the  two  scenes  Is  proportional  to 
the  correlation.  When  the  two  are  exactly  aligned,  the  overlap  (correlation) 
is  a  maximum  and  when  they  are  completely  separated,  the  overlap  (correlation) 
is  a  minimum.  If  light  is  shone  through  these  transparencies,  the  transmitted 
intensity  is  a  maximum  when  the  scenes  exactly  align.  This  reasoning  carries 
over  into  the  Fourier  transform  plane  in  much  the  same  way. 


D.  Sensitivity  to  Rotation 

As  a  simple  example,  consider  a  Ronchi  ruling  as  the  input  scene  to 
the  optical  correlator.  A  Ronchi  ruling  in  its  usual  capacity  is  a  series  of 
opaque  straight  parallel  lines  on  a  transparent  substrate;  the  width  of  the 
opaque  lines  being  equal  to  the  separation  between  them.  The  resulting 
matched  filter  consists  of  a  series  of  exposed  spots  -  roughly  circular  -  on 
either  side  of  a  strong  zeroth  order,  very  similar  to  what  would  be  expected 
in  the  far  field  of  a  diffraction  grating  illuminated  by  a  coherent  light.  A 
more  exact  description  is  given  in  Appendix  A.  If  the  exposed  plate  looks 
similar  to  Figure  13,  the  distance  from  the  center  of  the  zeroth  order  to  the 
center  of  either  of  the  closest  dark  circles  is  given  by 


xf  =  nXf  , 


(63) 


where  P  is  the  fundamental  spatial  frequency  of  the  Ronchi  ruling,  equal  to 
l/(2a),  where  a  is  the  spacing,  X  is  the  wavelength  of  the  illumination  and  f  is 
the  focal  length  of  the  Fourier-transforming  lens.  If  the  matched  filter  just 
made  is  addressed  with  the  same  Ronchi  ruling  input,  it  is  desired  to  deter¬ 
mine  how  the  correlation  will  change  when  the  input  is  rotated.  Figure  14 
illustates  this  phenomena.  The  correlation  should  be  dependent  on  the  area  of 
overlap  of  the  two  circles  in  the  Fourier  transform  plane.  This  area  is  given 
by :  [41] 


n2 

-1 

/  d  \ 

d  \ 

R2  d2 

R 

COS 

(  2R  ) 

2  ’ 

R  -  X  _ 

(64) 

Thus,  writing  d«  Xf0,  the  correlation  should  be  proportional  to:  [42] 


C 


(65) 


where  rtR^  is  a  normalization  factor.  A  plot  of  this  equation  Is  given  in 
Figure  15.  This  may  seem  like  a  trivial  example  but  it  contains  information 
about  how  the  correlation  intensity  changes  with  the  rotation  angle  and  gives 
an  indication  of  the  role  played  by  the  spatial  frequency  content  of  the 
scene.  The  variable  Xf  in  Equation  (65)  is  related  to  the  spatial  frequency: 


o  = 


'  f 


(66) 


Thus,  Equation  (65)  may  be  rewritten  as: 


f  = 

28.26  cm 

ft  = 

11.7  /  cm 

R  = 

. 007  cm 

x  = 

6.328  x  10 

xf= 

.021  cm 

.54  Radians 

31  Degrees 

Figure  15. 
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A  plot  of  equation  65  —  predicted  correlation 
versus  rotation  of  input  scene. 
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This  is  illustrated  in  Figure  16  for  several  different  values  of  P.  As  one 
can  see,  the  width  of  the  correlation  curve  broadens  when  the  spatial  fre¬ 
quency  decreases.  This  will  prove  useful  in  understanding  some  experimental 
results  given  later. 


E.  Sensitivity  to  Scale 

In  the  past,  one  of  the  problems  with  using  Fourier  transform 
matched  filters  has  been  the  sensitivity  of  the  correlation  to  a  scale  change. 
This  scale  change  might  be  brought  about  by  an  airplane  changing  altitude  or 
the  target  moving  in  the  field  of  view  toward  or  away  from  the  detector 
device.  Other  researchers  have  observed  a  50  percent  loss  in  correlation  for 
a  scale  change  of  as  little  as  0.5  percent  [42].  A  simple  model  of  this 
effect  is  possible.  Consider  an  input  scene  having  a  single  spatial  frequency 
content  P,  then  the  location  in  the  Fourier  transform  plane  of  the  first  order 
exposed  spot  will  be  a  distance: 


Xf  =  PAf  (68) 

from  the  zeroth  order.  See  also  the  photograph  of  Figure  29  and  the  sub¬ 
sequent  calculations.  A  change  in  scale  of  the  input  scene  produces  an  effec¬ 
tive  linear  change  in  the  spatial  frequency  content  of  the  scene  now 
addressing  the  matched  filter.  This  produces  a  linear  shift  in  the  location 
of  the  Fourier  transform.  Let  the  new  distance  from  the  zeroth  order  be  given 
by: 

=  q'\ f  (69) 


and  let 


op'  =  P 


(70) 


where  a  is  the  scale  change  (magnification)  factor.  If  we  assume  the  Fourier 
transform  elements  to  be  roughly  circular  -  see  Figure  17,  then  the  area  of 
overlap  should  be  proportional  to  the  correlation,  C: 


A 


o 


2 


(71) 


Thus,  the  correlation  as  a  function  of  magnification  should  be: 


(72) 
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where  i&2  is  a  normalization  factor.  This  is  illustrated  in  Figure  18.  Note 
that  a  magnification  of  about  1.37  produces  a  50  percent  loss  in  correlation. 
This,  of  course,  depends  on  the  spatial  frequency  content  much  the  same  as  the 
curves  in  Figure  16. 


F.  Sensitivity  to  Tilt 


Using  the  previous  analysis  it  will  be  fairly  simple  to  analyze  the 
effect  of  tilting  the  input  scene.  In  reality  this  might  occur  any  time  an 
off-normal  approach  would  be  preferred  or  if  the  object  began  to  turn  on  a 
longitudinal  axis.  Tilting  the  input  scene  (or  the  correlator)  has  the  effect 
of  changing  the  spatial  frequency  content  in  one  direction.  Assume  the  input 
scene  is  a  Ronchi  ruling  having  a  spacing  of  a  .  Tilting  the  ruling  with  the 
lines  as  an  axis  reduces  the  effective  width  by  a  factor  of  cos  <J>-  The  fun¬ 
damental  spatial  frequency  of  the  ruling  is  given  by 


Q  =  — 
2a 


(73) 


The  effective  spatial  frequency  of  the  tilted  ruling  is: 


O  3 


2a  cos 


(74) 


Thus 


_1_ 

2a 


Xf  =  —  \f 


(75) 


and 


Xr  =  —  sec  <•> 
f  2a 


(76) 


Substituting  this  into  Equation  (71)  gives: 


A  plot  of  this  is  given  in  Figure  19.  Notice  that  for  this  simple  model,  a  38 
degree  tilt  is  allowed  before  the  correlation  drops  by  50  percent.  The  corre¬ 
lation  does  not  seem  to  be  nearly  as  sensitive  to  tilt  as  it  is  to  rotation 
and  scale  change  (magnification). 

The  three  previous  models  of  rotation,  scale,  and  tilt  are  given  as 
illustrations  only.  In  reality  an  input  scene  contains  many  spatial  frequen¬ 
cies  and  usually  there  is  not  a  single  largely-dominate  one.  A  direct  com¬ 
parison  between  these  models  and  experimental  data  is  not  useful  except  for 
highly  specialized  scenes  which  rarely  occur  in  the  real  world.  These  models 
do,  however,  illustrate  what  might  be  expected  in  general. 
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G.  The  Multi-Foci  Hololens 

The  details  of  the  production  of  the  multi-foci  hololens  are  given  in 
papers  by  Liu  and  Duthie  [33]  and  Gregory  and  Liu  [34]  and  are  not  repeated 
here.  A  brief  description  is  given  which  highlights  the  analysis  and  experi¬ 
mental  procedure.  The  essential  element  is  a  screen,  patented  by  Liu,  which 
will  produce  a  fairly  equal  Intensity  5X5  array  of  spots  each  of  which  contains 
the  information  of  the  incident  beam.  This  screen  was  developed  after  much 
experimentation  and  many  computer  simulations.  A  photomicrograph  of  a  few  unit 
cells  of  this  screen  is  given  in  Figure  20.  Crudely  speaking,  the  screen  con¬ 
sists  of  an  array  of  crossed  broken  Ronchi  rulings.  With  the  aid  of  carefully 
placed  neutral  density  filters,  this  screen  will  produce  a  5X5  array  of  spots 
which  serves  as  the  input  to  the  optical  system  given  in  Figure  21.  The  inci¬ 
dent  beam  I  (a  collimated  argon  ion  laser)  impinges  upon  the  contact  screen  S 
discussed  above;  (the  neutral  density  filters  are  not  shown),  thus,  producing 
an  array  of  5X5  bright  spots  as  well  as  other  weaker  orders.  A  lens  Lj  is 
placed  a  focal  length  f^  from  the  screen  S.  In  the  focal  plane  of  this  lens  a 
square  aperture  A  is  placed  to  block  all  but  the  5X5  matrix.  This  is  then  the 
input  to  a  standard  holographic  arrangement.  The  placement  of  lens  L2 
determines  the  magnification  of  the  resulting  holographic  element.  At  a 
distance  f2  (the  focal  length)  from  the  lens,  a  dichromated  gelatin  pho¬ 
tographic  plate  is  placed.  A  reference  beam,  derived  from  the  original  argon 
ion  beam,  is  allowed  to  Interfere  with  the  object  beam  at  the  holographic  plate 
F.  The  exposed  plate  is  then  developed  and  sealed  with  a  thin  glass  cover 
plate  since  dichromated  gelatin  holograms  are  known  to  be  humidity-sensitive 
[43].  The  holographic  plate,  when  reilluminated,  will  thus  produce  a  5X5  array 
of  Fourier-transform  elements  off-axis  and  along  the  path  of  the  original 
reference  beam;  see  Figure  22.  Figure  23  is  an  actual  photograph  of  the  array  - 
magnified  about  five  times. 

The  array  is  produced  in  focus  a  distance  f^  away  from  the  plate.  If 
the  distance  dQ  in  Figure  21  is  called  the  object  distance  then  the  image 
distance  is  given  by: 

di  =  f2  +  fh  (78) 

so  that  the  geometrical  thin  lens  equation  is: 

_L  .  JL  .  _L  •  (79) 

d0  dt  f2 

Thus,  the  focal  length  of  the  hololens  can  be  written  as: 


_1_ 


The  magnification  is  given  by: 


M  = 


d 
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f  2  +  fh 


(80) 
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Therefore,  It  Is  possible  to  tailor  a  hololens  to  meet  the  specific  require¬ 
ments  of  the  experiment.  The  hololens  used  in  this  investigation  had  a  focal 
length  of  28.26  centimeters  and  produced  an  off-axis  array  of  5X5  spots  with  a 
separation  of  2.7  mm  between  them.  This  combination  proved  to  be  quite 
satisfactory. 


III.  EXPERIMENTAL  PROCEDURES 


A.  The  Experimental  Arrangement  and  Operation 

The  basic  experimental  arrangement  used  to  make  matched  filters  has 
not  changed  appreciabiy  in  20  years.  The  Vander  Lugt  method,  as  outlined  in 
Section  I,  is  still  the  most  widely  used  in  optical  data  processing.  The  fun¬ 
damental  elements  are  a  Fourier-transforming  lens  (or  mirror),  a  coherent  light 
source  and  a  medium,  such  as  a  photographic  plate,  for  recording  the  filter. 

The  input  scene  may  be  a  transparency  or,  in  this  case,  a  coherent  image  pro¬ 
duced  by  a  LCLV. 

Figure  24  gives  the  experimental  arrangement  as  set  up  on  a  granite 
table.  The  image  on  TV  monitor  1  is  transmitted  from  TV  camera  1  which  views 
the  scene  mounted  on  a  table  which  has  rotation  and  tilt  capabilities.  The 
image  produced  by  TV  screen  1  is  focused  onto  the  "write"  side  of  the  Hughes 
LCLV.  The  illumination  of  the  "read"  side  of  the  LCLV  is  provided  by  a  15  MW 
HeNe  laser  ( Spectraphysics  model  Stabilite).  After  going  through  a  beam¬ 
expanding  and  collimating  arrangement,  the  beam  is  split  using  a  dielectric 
coated  beamsplitter  into  two  reasonably  equal  intensity  beams.  One  of  these 
beams  becomes  the  reference  beam;  the  other  strikes  the  LCLV  causing  the 
polarization  of  the  beam  to  be  spatially  modulated  according  to  the  intensity 
modulation  provided  by  the  television  screen  input.  The  details  of  the  opera¬ 
tion  of  the  LCLV  are  given  in  a  later  section. 

The  beam  containing  spatial  potarlzation  information  about  the  input 
scene  travels  back  from  the  LCLV  through  the  beamsplitter  and  on  to  polarizer 
1.  The  function  of  this  polarizer  is  to  turn  the  polarization-modulated  beam 
back  into  an  intensity-modulated  beam.  The  purpose  of  the  apparatus  described 
to  this  point  is  to  produce  a  coherent  image  of  the  input  scene.  The  holoiens 
performs  two  functions:  the  coherent  input  image  is  split  into  25  identicai 
images  and  each  of  these  images  is  Fourier-transformed.  The  details  of  this 
holographic  element  were  discussed  in  Section  II.  A  photographic  plate  (Kodak 
649F  or  120)  is  placed  a  focal  length  (28.26  cm)  behind  the  lens.  An  exposure 
is  made  by  removing  the  beam  block  from  the  reference  beam  and  activating  the 
shutter  for  a  specified  amount  of  time.  The  purpose  of  polarizer  2  is  to 
Insure  that  the  object  and  reference  beam  are  of  the  same  linear  polarization. 

The  mask  shown  In  Figure  24  allows  only  one  of  the  elements  of  the 
array  to  pass  along  with  the  reference  beam.  The  location  of  the  mask  is 
computer-controlled  for  accurate  positioning.  After  exposing  as  many  of  the 
elements  as  desired,  the  plate  is  removed  and  developed  using  standard  tech¬ 
niques  [44].  The  developing  procedure  is  given  in  Appendix  B.  The  developed 
plate  is  then  placed  back  In  the  system  In  Its  original  position.  The 
reference  beam  is  then  blocked  and  the  mask  removed  (for  parallel-addressing). 
The  object  beams  fail  on  the  matched  filters  just  created  and  if  any  one  of  the 
filters  correlates  with  the  input  scene  now  on  the  TV  screen,  a  reconstructed 
reference  beam  wiii  result  from  the  plate.  This  beam  is  Fourier-transformed  by 
lens  L2  onto  an  RCA  CCD  camera  (Model  TC  1160)  and  viewed  as  a  localized  bright 
spot  on  TV  monitor  2.  The  Colorado  video  analyzer  (Model  321)  translates  the 
brightness  of  the  spot  into  a  voltage  which  is  input  to  a  chart  recorder 
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(Hewlett  Packard  Model  680).  The  correlation  intensity  may  also  be  displayed  on 
an  oscilloscope  directly  and  photographed  if  desired.  The  input  scene  can  then 
be  rotated,  tilted,  translated,  or  changed  completely  as  the  chart  recorder 
displays  the  change  in  intensity  of  the  spot. 

B.  The  Liquid  Crystal  Light  Valve  (LCLV) 

Figure  25  gives  a  schematic  drawing  of  the  LCLV  used  in  this 
research.  Details  of  the  fabrication  and  testing  of  this  type  device  are  given 
in  several  papers  [45]  [46).  The  LCLV  used  in  this  research  was  made  by  the 
Hughes  Corporation.  In  operation  [28),  incoherent  light  from  a  television 
monitor  or  other  source  passes  through  the  transparent  glass  face  plate  and 
through  the  transparent  conductive  electrode  (a  tin  oxide  indium  oxide  glass) 
and  strikes  the  photosensor  composed  of  a  layer  of  cadmium  sulfide.  This 
material  becomes  conductive  when  the  light  areas  of  the  image  are  present.  The 
charges  supplied  by  the  power  source  collect  on  the  surface  of  the  dielectric 
mirror  corresponding  to  the  conductive  areas  of  the  photosensor,  thus,  creating 
an  electric  field  across  the  twisted  nematic  liquid  crystal.  This  liquid 
crystal  has  been  preoriented  using  the  liquid  crystal  alignment  films  which 
consist  of  linearly  scratched  surfaces  so  that  the  nematic  liquid  crystal 
"needles"  go  through  a  45  degree  twist  from  one  surface  to  the  other.  The  pref¬ 
erentially  applied  electric  field  causes  selected  areas  of  the  liquid  crystal 
to  untwist  proportional  to  the  strength  of  the  localized  field  and  try  to  align 
with  the  applied  field.  The  plane-polarized  coherent  projection  light  enters 
through  a  transparent  glass  substrate,  passes  through  the  liquid  crystal  and  is 
reflected  back  by  the  dielectric  mirror.  Selected  areas  of  the  beam  will 
undergo  a  polarizatrion  change  provided  by  the  preferential  twist  and  the 
intrinsic  birefringence  of  the  liquid  crystal:  Dark  areas  of  the  incoherent 
"write"  image  produce  areas  in  the  "read”  beam  that  are  plane-polarized.  Light 
areas  produce  (in  general)  elliptically-polarized  areas  in  the  "read"  beam. 

Thus,  the  reflected  beam  is  polarization-encoded  with  the  original  incoherent 
image.  The  resulting  coherent  image  may  then  be  decoded  using  a  linear 
polarizer,  producing  an  image  with  a  continuous  intensity  variation.  Complete 
details  of  the  operation  of  this  device  are  given  by  Grinberg  [27], 

Since  its  production,  the  Hughes  LCLV  has  been  tested  and  analyzed  by 
several  authors  [47]  [48],  Resolution  has  been  one  of  the  main  problems  along 
with  response  time.  The  resolution  of  the  LCLV  used  in  this  investigation  was 
measured  using  a  series  of  Ronchi  rulings  consisting  of  flat  black  lines  on  a 
white  background.  Tt  was  found  that  slightly  less  than  20  lines  per  millimeter 
(lines/mm)  could  be  distinguished.  This  was  partly  due  to  the  resolution  of 
the  television  monitor  used  as  the  input.  Assuming  512  lines  across  a  10-inch 
television  screen  focused  to  a  1-inch  image  gives  a  needed  resolution  of  20 
lines/mm.  Hughes  quotes  the  LCLV  as  having  a  resolution  of  60  lines/mm  and  a 
modulation  transfer  function  of  50  percent  at  this  frequency.  If  this  is  true, 
it  seems  the  television  monitor  is  the  limiting  piece  of  equipment. 

The  reaction  time  was  not  measured  in  this  investigation  but  is 
quoted  by  Hughes  as  being  25  ms.  Other  characteristics  of  the  LCLV  are  given 
in  Table  1  [28].  Even  though  the  LCLV  certainly  has  its  limitations,  it  is 
still  the  overall  best  incoherent-to-coherent  converter  available  today. 

Research  is  underway  at  the  Army  Missile  Laboratory  in  an  attempt  to  improve 
its  performance. 
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Figure  25.  The  Liquid  Crystal  Light  Valve  (LCLV) 


TABLE  1.  Characteristics  of  the  Hughes  Liquid  Crystal  Light  Valve  (LCLV) 


Aperture  Size 

1  in. 2  (l  in.  X  1  in.) 

Sensitivity 

160  nw/cm2  at  5250A 

Resolution 

60  lines/mm  at  50%  MTF 

Contrast 

>  100:1 

Grayscales 

9 

Speed 

Excitation 

10  msec 

Extinction 

15  msec 

Reflectivity 

>  90% 

Optical  Quality 

<  X/4  at  6328A 

Voltage 

6V(rms)  @  10  KHz  (variable) 

Other  light  modulators  include  the  Pockel's  Read-Out  optical  Modu- 
ator  (PROM)  and  various  photoplastic  devices.  PROM  devices  are  fabricated  from 
several  electro-optical  crystal  materials  such  as  ZnS,  ZnSe,  and  Bi^SiC^o  and 

convert  an  incident  intensity  pattern  to  a  voltage  pattern.  The  areas  of  the 
crystal  not  experiencing  a  voltage  act  like  a  half-wave  retardation  plate  to 
the  coherent  "read”  light.  Thus,  the  coherent  image  is  polarization-encoded 
and  an  analyzer  (polarizer)  can  convert  it  back  into  an  intesity-encoded  image. 
This  device  shows  great  promise  but  at  the  present  time  is  quite  expensive  and 
cannot  maintain  the  coherent  image  for  very  long  when  even  moderately  intense 
illumination  is  used.  Photoplastics  are  phase  recording  materials  with  the 
input  intensity  being  recorded  as  deformations  on  a  pliable  surface.  At  pres¬ 
ent,  these  devices  are  quite  slow  and  not  very  sensitive  [39]  (p.  133). 

C.  Photographic  Techniques 

Photographic  techniques  are  an  integral  part  of  producing  holographic 
matched  filters.  In  this  investigation  Kodak  high  resolution  holographic 
plates  (649F  or  120)  were  used  exclusively.  These  plates  are  quoted  as  having 
a  resolution  of  better  than  2000  lines/mm  [49].  The  high  resolution  film 
allows  separation  between  spatial  frequencies  that  are  very  similar. 

Figure  26  is  a  photograph  of  a  magnified  matched  filter.  The  input 
scene  was  an  aerial  photograph  of  Huntsville,  AL  discussed  later  in  Section  IV. 
The  magnification  in  this  photograph  is  about  100;  thus,  the  actual  diameter  of 
the  filter  is  roughly  .035  cm.  The  minimum  sensitivity  of  Kodak  649F  is  about 
8  X  10“5  joules  per  square  centimeter  (J/cm^)  and  about  4  X  10“^  J/cm^  for  type 
120  [49].  The  exposed  area  is  roughly  1  X  10“^  cm2  (assuming  a  circle  of 
diameter  .035  cm).  The  average  power  of  each  of  the  elements  was  typically 
.15  \U. 


Figure  26.  A  photograph  of  a  matched  filter  of  an 
aerial  map  of  Huntsville,  AL  (M  =  100) 
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(83) 


Thus,  the  exposure  time  needed  to  reach  minimum  sensitivity  is  given  by: 
train  =  ~~A 

where  tmin  is  the  minimum  exposure  time,  P  is  the  power,  p  the  minimum 
sensitivity,  and  A  is  the  area  being  exposed.  Using  the  previously  given  rough 
values,  tmin  is  found  to  be  about  2  seconds  for  649F  film  and  about  4  seconds 
for  type  120.  In  practice,  exposure  times  of  10  to  12  seconds  were  used  to 
obtain  the  best  matched  filters.  The  calculated  times  served  as  a  good 
starting  point  to  begin  the  investigation.  The  emulsion  thickness  of  each  of 
the  films  was  about  17  |jm  for  649F  and  about  6  pm  for  type  120  [50].  Six  pm, 
is  considered  an  average  thickness  for  photographic  emulsions.  The  film  was 
developed  using  well  established  procedures:  6  minutes  in  slightly  agitated 
D-19  developer,  3  minutes  in  Kodak  stop  bath,  and  1  minute  in  Hunts  Aerofix 
fixer.  The  plate  was  then  placed  in  a  water  bath  for  at  least  5  minutes, 
rinsed  with  absolute  ethanol  and  allowed  to  air-dry.  Complete  details  are 
given  in  Appendix  B. 

D.  Multiply  Exposed  Filters 

Previously  it  was  mentioned  that  Christensen  was  able  to  obtain  good 
correlation  of  eight  targets  addressed  in  series  by  multiply  exposing  the  pho¬ 
tographic  plate  [31].  This  was  also  used  in  this  investigation;  results  are 
given  in  Section  IV.  A  phenomena  called  the  intermittency  effect  should  be  men¬ 
tioned  at  this  point:  several  exposures  are  not,  in  general,  equal  to  a  single 
exposure  of  the  same  total  duration  [51].  A  paper  by  Johnson  [52],  dealing 
with  multiply  exposed  holograms  gave  data  that  indicated  that  the  exposure 
times  used  should  be  roughly  a  function  of  the  square  of  the  number  of  expo¬ 
sures  made: 

eN  =  AN2  +  B  (84) 

where  A  and  B  are  constants  and  N  is  the  order  of  the  exposure.  ejj  is  the 
exposure  time  for  the  exposure.  Biedermann  [53]  suggests  the  exposure  time 
be  increased  about  12  percent  over  the  previous  one  for  each  successive  expo¬ 
sure  using  Scientia  10E70  plates  which  are  similar  to  Kodak  649F.  In  the  pre¬ 
sent  investigation,  multiple  exposures  were  made  by  rotating  the  input  scene  by 
computer  control  between  exposures.  The  exposure  times  required  were  also 
programmed  into  the  computer  as  necessary. 

E.  Equipment 

The  following  is  a  list  of  the  equipment  used  throughout  this  investi¬ 
gation  along  with  appropriate  identification  and  serial  numbers  where 
applicable.  The  USA  designation  indicates  a  local  property  control  number: 

1.  Helium-Neon  laser,  Spectraphysics  Model,  Stabilite,  S/N  CJ-111 

(USA). 

2.  Liquid  crystal  light  valve,  Hughes  Research,  45°  biphenyl 
nematic,  S/N  10962,  CJ  137  (USA). 


3.  Television  camera,  CCD,  RCA  Electronics,  Model  TC-1160,  S/N 

B001669. 

4.  Television  camera,  vidicon,  RCA  Electronics,  Model  TC-1005,  S/N 
CJ  043  (USA). 

5.  Image  digitizer,  Colorado  Video  Model  321,  S/N  231. 

6.  Hololens,  Lumin  Electro-Optics  (H.  K.  Liu),  Model  LUMIN025. 

7.  Television  monitor,  Videotek,  Model  VM12PR,  S/N  1282118B. 

8.  Television  monitor,  Sony  Model  CVM-194,  S/N  015612. 

9.  Computer,  Hazeltine,  Model  1500,  S/N  CJ266  (USA). 

10.  Translation  stages,  Ealing  (2  each),  Model  352427  and  352419,  S/N 
CJ  273  (USA)  and  CJ  103  (USA). 

11.  Ealing  rotation  stage,  Model  344614,  S/N  CJ  141  (USA). 

12.  Isolation  table.  Granite  Do-All,  S/N  2486-3. 

13.  Power  meter,  NRC,  Model  880,  S/N  1010178. 

14.  Photographic  plate,  Kodak  649F  and  type  120. 

15.  Various  lenses,  polarizers  and  beamsplitters  supplied  by  several 

vendors. 

16.  Power  supply,  Hewlett-Packard,  Model  8116A,  S/N  CJ  089  (USA). 

17.  Rotation  and  tilt  stage.  Scientific  Atlanta,  Model  4116-10,  S/N 
CJ  051  (USA). 

18.  Oscilloscope,  Tektronix,  Model  7603,  S/N  CJ  162  (USA). 

19.  Chart  recorder,  Hewlett-Packard  Model  680 

20.  Digital  voltmeter,  Keithley,  Model  179. 
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EXPERIMENTAL  RESULTS 


A.  Testing  of  Hoioiens 

The  hoioiens  used  in  this  investigation  was  made  by  H.  K.  Liu  of  the 
University  of  Alabama,  Tuscaloosa.  A  lens  of  this  type  had  never  been  used  or 
tested  before  so  two  basic  experiments  were  devised  to  measure  its  performance. 
The  screen  used  by  Liu  and  Duthie  had  several  problems  associated  with  it,  the 
major  one  being  the  lack  of  uniformity  in  intensity  of  the  array  produced  by 
the  element  [33].  The  first  test  was  to  photograph  the  array  using  a  long 
exposure  and  look  for  extraneous  diffraction  orders.  The  photograph  is  shown 
in  Figure  27.  One  can  see  that  the  extra  orders  are  visible  but  are  down  in 
intensity  several  orders  of  magnitude  when  compared  with  the  25  primary  ele¬ 
ments  of  the  array;  thus,  the  "cross-talk”  was  not  expected  to  be  a  problem. 

The  separation  between  the  elements  was  found  to  be  a  regular  2.7  mm.  The 
focal  length  of  the  hoioiens  was  measured  to  be  28.26  cm.  In  order  to  have 
"cross  talk"  between  elements  It  would  be  necessary  to  have  a  scene  with  a  spa¬ 
tial  frequency  component  which  would  extend  half  the  distance  between  elements, 
thus,  giving  overlapping  exposures.  Equations  (66)  predicts  this  spatial  fre¬ 
quency  to  about  75/cm  which  is  almost  four  times  the  highest  fundamental  spa¬ 
tial  frequency  that  can  be  obtained  using  a  10-inch,  512  line  standard 
television. 


As  previously  mentioned,  the  Intensity  distribution  in  earlier 
diffracting  holographic  elements  was  quite  nonuniform.  The  power  contained  in 
each  element  of  the  present  array  was  carefully  measured  using  a  detector  (NRC 
Model  880)  tuned  to  the  HeNe  wavelength.  The  results  of  this  measurement  are 
given  in  Table  2.  As  one  can  see,  there  is  still  a  variation  in  intensity  over 
the  field,  but  the  uniformity  is  greatly  improved  over  that  measured  by  Liu  and 
Duthie. 


TABLE  2.  Power  Distribution  in  the  Image  Plane  of  the  Hoioiens. 


Note:  Power  distribution  (in  pw)  of  each  of  the  25  focuses  of  the 
holographic  element  with  a  uniform  LCLV  input.  LCLV  is  driven  at  f  *  1400  Hz 
and  V  ■  6.5  volts. 


The  voltage  and  current  used  to  drive  the  LCLV  Is  typical  of  that  used 
throughout  this  investigation.  It  is  Interesting  to  observe  that  the  average 
power  over  all  25  elements  is  .54  yW  which  is  the  value  of  the  power  measured 
at  location  (-1,2).  This  experiment  was  repeated  several  times  with  different 
adjustments  of  the  LCLV  with  the  same  result;  position  (-1,2)  was  the  average 
power  of  the  array.  The  standard  deviation  c  was  calculated  to  be  .18  pW  for 
the  25  data  elements  of  Table  2.  The  variation  in  intensity  of  the  elements 
generally  caused  no  major  problems  until  the  last  phase  of  the  research  when 
many  filters  were  stored  at  each  of  the  elements.  The  exposure  times  and  beam 
ratios  became  critical  at  this  point  and  a  very  uniform  array  would  have  made 
for  better  correlations. 

The  efficiency  of  the  hololens  was  next  investigated.  It  is  known 
that  holograms  stored  using  dichromated  gelatin  plates  can  have  efficiencies  of 
70  percent  or  better  [54],  The  hololens  made  by  Liu  and  used  in  this  investi¬ 
gation  was  quite  disappointing  in  this  respect.  An  incident  power  of  104  pW  in 
a  well-collimated  HeNe  beam  resulted  in  a  total  power  of  5  pW  in  the  5X5  array. 
This  is  an  efficiency  of  about  5  percent.  The  remaining  power  went  mostly  into 
the  zeroth  (undiffracted)  order.  The  laser  used  for  this  entire  work  was  a  15 
mW  Spectraphysics  (Stabilite  Model).  Thus,  even  with  a  low  efficiency  hololens 
it  was  possible  to  obtain  enough  power  to  produce  the  exposures  needed. 

Because  of  the  low  efficiency,  however,  some  exposure  time  were  inordinately 
long  (15  seconds)  which  can  certainly  cause  problems  if  vibrations  of  the 
system  are  present.  This  problem  was  not  observed,  however. 

The  next  test  was  to  see  if  the  screen  actually  performed  that  task 
of  a  lens  under  the  same  circumstances.  That  is,  does  each  spot  represent  the 
Fourier  transform  of  the  coherent  input?  According  to  Goodman  [38]  (p.  86), 
the  distribution  of  the  field  in  the  back  focal  plane  for  a  transparency 
(aperture,  etc.)  located  in  the  front  focal  plane  illuminated  with  a  coherent 
plane  wave  is: 


U(xf ,Yf) 


—  On 


(85) 


where  A  is  the  amplitude  of  the  Incident  plane  wave,  f  is  the  focal  length  of 
the  lens  and  \  is  the  wavelength  of  the  light.  (x0*yo)  represents  the  front 
focal  plane  while  (xf.y^)  denotes  the  Fourier  transform  plane  and  t0  (x0,y0)  is 
the  transmission  function. 

For  a  circularly  symmetric  aperture  a  transformation  of  coordinates 
is  possible,  resulting  in  [38]  (p.  16) 
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where 


P  = 


_1_ 

Xf 


2 

yf 


(87) 


is  the  radial  spatial  frequency,  is  a  Bessel  function  of  the  first  kind, 
order  one  and  R  is  the  radius  of  the  circular  aperture.  The  intensity  is  then 
proportional  to  U(  p)  times  its  complex  conjugate. 


I(p) 


A2tt2R4 

X2f2 
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(88) 


The  first  zero  in  occurs  when 
2  pR  =  1.22 


(89) 


or 


2R 

Xf 


1.22  . 


(90) 


Looking  only  at  the  Xf  dependence  because  of  circular  symmetry: 


2R 

\  f  Xf  =*  1.22 

or,  solving  for  xf 


(91) 


xf 


1.22  -f 
D 


where  D  =  2R  is  the  diameter  of  the  aperture. 


(92) 


If  the  hololens  actually  behaves  as  well  as  a  lens,  then  this  equation 
should  be  valid  when  a  circular  aperture  is  placed  a  focal  length  before  the 
element  (28.26  cm).  Using  an  aperture  of  about  2.5  mm  (an  adjustable  iris)  and 
a  wavelength  of  .6328  pm,  the  pattern  of  Figure  28  results.  This  is  a 
photograph  of  the  exposed  and  developed  Kodak  649F  plate  used  to  record  the 
intensity  pattern.  The  distance  to  the  center  of  the  first  dark  ring  is  pre¬ 
dicted  to  be 


xf 


1.22  (.6328  X  10~6M)  (.2826  M) 
.0025  M 


(93) 


Xf  =  .0089  cm.  (94) 

A  measurement  yielded  .0080  cm  which  is  reasonable  agreement  (10  percent  error) 
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A  Ronchi  ruling  type  input  scene  was  then  used  as  the  input  to  the 
television  camera  of  Figure  24.  The  ruling  consisted  of  equally  spaced  flat 
black  lines  on  a  white  background.  The  lines  were  1/4-inches  wide  and 
1/4-inches  apart.  A  photograph  of  the  Fourier  transform  produced  is  given  in 
Figure  29.  The  Fourier  transform,  in  this  case,  was  viewed  via  a  television 
camera,  with  associated  optics,  placed  in  the  focal  plane.  The  image  was  then 
viewed  on  a  television  monitor  and  photographed.  The  separation  between  spots 
should  be  (from  Section  II): 


or 


xf  -  QXf  (96) 

where  Q  is  the  spatial  frequency,  X  is  the  wavelength  (HeNe)  and  f  is  the  focal 
length  of  the  hololens  (28.26  cm).  The  fundamental  spatial  frequency  of  the 
input  scene  is  2  per  inch  or  .79  per  cm.  In  this  case  a  magnification  factor  of 
1/8  is  involved  in  transferring  the  image  from  the  TV  camera  to  the  output  of 
the  LCLV  so  that  the  actual  spatial  frequency  is  about  6.3  per  cm.  The  above 
equation  thus  predicts  the  separation  between  the  zeroth  and  first  order  to  be: 

xf  =  (6.3  cm~l ) (6.328  X  10"5  cm)  (28.26  cm).  (97) 

Thus,  Xf  should  be  about  .011  cm.  The  magnification  of  the  photo  in  Figure  29 
was  found  to  be  about  80.  This  gives  an  average  measured  separation  of  about 
.012  cm,  which  is  within  10  percent  of  the  calculated  value. 

The  pattern  of  Figure  29  is  not  very  distinct  due  to  the  low  spatial 
frequency  of  the  input  scene.  It  might  be  compared  to  results  obained  when 
looking  at  the  diffraction  pattern  produced  by  a  few  slits.  The  greater  the 
number  of  slits,  the  more  distinct  the  pattern.  The  tests  thus  completed  showed 
that  the  hololens  does,  with  reasonable  accuracy,  perform  at  least  some  of  the 
tasks  of  an  ordinary  lens.  Difficulties  did  show  up  later  and  will  be  discussed 
in  this  Section. 

B.  The  K  Ratio 

In  Section  III  the  experimental  procedure  outlining  the  making  of  a 
matched  filter  was  given.  A  detail  not  mentioned  was  that  of  determining  the 
optimum  K  ratio.  The  K  ratio  is  defined  as  the  ratio  of  the  reference  beam 
intensity  to  the  object  beam  intensity  [25].  It  might  seem  obvious  that  this 
ratio  should  be  one  -  the  reference  beam  intensity  being  equal  to  the  object 
beam  intensity.  It  seems  that  this  ratio  should  optimize  the  interference  be¬ 
tween  the  two  beams  at  the  photographic  plate.  However,  the  intensity  at  the 
photographic  plate  corresponding  to  one  of  the  spatial  frequency  components  of 
the  input  scene  will  not,  in  general,  be  equal  to  the  intensity  contribution 
due  to  a  different  spatial  frequency  component  of  the  input  scene.  Thus,  it  is 


impossible  to  determine  a  single  K  value  that  will  be  valid  for  ail  input  scenes 
This  value  must  be  determined  experimentally.  This  knowledge  will  allow  a 
matched  filter  to  be  tailored  to  a  specific  frequency  component  of  the  scene 
if  desired.  As  an  example,  if  only  the  outline  of  an  object  is  of  importance, 
such  as  an  airplane,  then  the  filter  can  be  made  using  a  K  ratio  tuned  to  maxi¬ 
mize  the  interference  in  the  low  spatial  frequency  range.  If,  however,  the 
problem  is  to  distinguish  between  a  plane  with  windows  and  a  plane  with  no  win¬ 
dows,  the  higher  spatial  frequencies  must  be  emphasized.  In  this  investigation 
a  scene  comtaining  both  high  and  low  spatial  frequency  components  was  used  to 
obtain  a  K  value  to  start  from  initially.  The  scene  was  an  aerial  photograph 
of  Huntsville,  AL,  shown  in  Figure  30.  As  one  can  see  there  are  both  fine  and 
coarse  details  present.  This  also  afforded  an  opportunity  to  study  the  expo¬ 
sure  times  that  would  be  needed.  Table  3  gives  the  results  for  a  series  of 
matched  filters  made  of  the  Huntsville  scene.  The  intensity  of  the  reference 
beam  was  changed  with  a  pair  of  polarizers. 


TABLE  3.  K-Ratio  Values. 
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t  (sec) 
Intensity 

(Arbitrary  Units) 

PR 

D 

9  sec 

12  s 

15  s 

18  s 

21  s 

0 

1.1 

1.34 

1.62 

1.47 

1.27 

0.98 

6  s 

9  s 

12  s 

15  s 

18  s 

2/1 

1.49 

1.56 

1.10 

0.60 

0.34 

3  s 

6  s 

9  s 

12  s 

15  s 

3/1 

1.12 

1.52 

0.83 

0.43 

0.20 

3  s 

6  s 

9  s 

12  s 

15  s 

5/1 

1.48 

0.40 

0.14 

0.00 

0.00 

3  s 

6  s 

9  s 

12  s 

15  s 

10/1 

0.62 

0.00 

0.00 

0.00 

0.00 

PQ  =  .15  pW 

NOTE:  The  intensities  of  the  correlation  output  vs  the  reference 

PR 

to  object  beam  ration  —  and  the  exposure  time  (seconds)  of 

*o 

a  single  matched  filter  of  the  photograph  of  the  city  of 


Huntsville. 


With  available  equipment,  it  was  only  possible  to  measure  power;  so 
in  order  to  get  an  estimate  of  the  intensity  of  the  beams,  a  mask  with  a  small 
( 1  mm)  hole  was  placed  over  the  detector  and  used  consistently  for  all 
measurements.  Thus,  the  value  of  K  in  this  work  may  not  be  compared  with  that 
given  by  other  researchers.  In  many  cases  authors  will  call  K  the  ratio  of  the 
total  power  in  the  reference  beam  to  the  total  power  in  the  Fourier-transformed 
object  beam.  This  K  will  be  very  large  on  occasions  when  the  reference  beam  is 
much  larger  than  the  object  beam  as  is  often  the  case.  This  method  of  deter¬ 
mining  K  did  not  seem  reasonable  for  this  investigation;  and  equipment  was  not 
readily  available  to  determine  K  as  Casasent  did  using  specialized  detectors 
and  micropositioning  techniques  [55],  thus,  a  compromise  was  sought.  Placing  a 
mask  over  the  detector  with  a  hole  in  it  large  enough  to  ensure  that  all  the 
spatial  frequency  components  would  pass  and  using  this  same  raask/detector 
arrangement  to  measure  the  power  in  the  reference  beam  seemed  to  be  a  reason¬ 
able  way  of  determining  a  "working  value"  of  K.  The  same  mask  and  detector 
were  used  throughout  this  investigation. 

The  intensity  of  the  object  beam,  PQ,  remained  constant,  as  it  was 
not  easily  changed  without  changing  the  quality  of  the  resulting  image.  The 
best  image  resulted  in  an  average  power  in  the  array  of  .  15pW  per  element.  The 
best  correlations  were  obtained  using  a  K  between  1  and  3  and  an  exposure  time 
of  about  12  seconds. 

C.  A  Single  Matched  Filter 

The  previous  tests  and  calculations  set  the  foundation  for  making  the 
first  matched  filters  themselves  using  the  hololens.  It  was  desired  to  test 
the  ability  of  the  hololens  against  a  glass  lens  in  making  matched  filters. 

The  hololens  was  replaced  by  a  high  quality  glass  lens  of  approximately  the 
same  focal  length  and  a  matched  filter  made  of  the  aerial  Huntsville  scene. 
Figure  32  gives  the  comparison  between  the  correlation  obtained  using  a  glass 
lens  and  that  obtained  with  one  element  of  the  hololens.  The  input  scene  was 
rotated  to  produce  the  curve  shown.  As  one  can  see,  the  hololens  performed 
well.  The  full  width  at  half  maximum  (FWHM)  is  about  2.8  degrees  for  the  holo¬ 
lens  and  slightly  less  for  the  glass  lens.  The  background  observed  using  the 
hololens  is  discussed  in  Paragraph  E  of  this  Section.  The  correlation  appeared 
as  a  very  bright  spot  on  the  television  monitor  and  was  easily  detected.  These 
particular  filters  were  made  using  a  K  value  near  one  and  an  exposure  time  of 
12  seconds  and  649F  film.  The  correlation  curve  is  quite  symmetric  about  zero 
degrees  rotation  -  the  position  of  the  scene  where  the  matched  filter  was  made. 
This  was  found  to  be  true  to  a  high  degree  in  all  cases  studied.  It  is  obvious 
that  the  signal-to-noise  ratio  is  enormous  for  this  single  matched  filter. 

This  is  not  the  case  however,  when  many  filters  are  stored.  This  is  demon¬ 
strated  in  the  following  sections.  The  matched  filters  themselves  were  very 
similar  to  the  photograph  of  Figure  26.  It  is  interesting  to  try  to  calculate 
the  highest  spatial  frequency  recorded.  In  Figure  26  the  most  distant  dark 
8 pot  is  at  a  distance  of  about  .055  cm.  Equation  (66)  predicts  this  to  corres¬ 
pond  to  a  spatial  frequency  of  about  30  per  cm.  This  is  somewhat  in  agreement 
with  Upatnieks  for  a  similar  input  scene  and  experimental  arrangement  [56]. 
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D.  Multiple  Exposures  -  Eight  Filters 


Other  researchers  have  shown  that  it  Is  possible,  using  a  conven¬ 
tional  glass  Fourier-transforming  lens,  to  store  several  matched  filters  at  a 
single  location  and  address  the  filter  in  series  [57].  A  test  was  thus  devised 
in  this  Investigation  to  examine  the  possibility  of  doing  the  same  thing  with 
the  hololens.  A  mask,  was  used  to  block  ail  but  one  element  of  the  5X5  array 
and  this  element  was  exposed  eight  times.  The  input  scene  was  again  the  aerial 
view  of  Huntsville.  In  this  instance,  2.5  seconds  was  used  for  each  exposure 
with  a  K  value  of  1.  The  resulting  correlations  are  given  in  Figure  33.  The 
intensities  are  weak  but  clearly  visible.  The  signai-to-noise  is  greater  than 
2  for  aii  filters.  All  eight  of  the  filters  have  about  the  same  FWHM  -  about  3 
degrees  of  rotation.  Very  little  hololens  background  noise  was  observed,  as  was 
the  usual  case  when  correlations  were  read  out  in  series.  Some  of  the  higher 
frequency  noise  superimposed  on  the  curves  was  electrical  in  nature.  The  sen¬ 
sitivity  of  the  equipment  had  to  be  increased  in  order  to  detect  the  weaker 
signals.  Some  of  the  noise  could,  no  doubt,  be  filtered  out  with  additional 
equipment.  This  was  not  attempted  in  this  investigation.  This  part  of  the 
research  did  show  up  to  eight  filters  can  be  stored  using  one  element  of  the 
hololens.  This  is  about  the  best  that  has  been  done  using  a  conventional  lens. 

E.  Twenty-Five  Filters  of  One  Scene  -  Hololens  Noise 

Early  in  the  Investigation  aii  25  elements  were  exposed  at  the  same 
time  with  a  single  input  scene.  This  afforded  the  first  observation  of  a 
problem  with  the  hololens.  Figure  34  shows  a  typical  correlation  curve  for  the 
rotated  aerial  view  of  Huntsville.  Notice  the  nonzero  baseline  even  after  a 
large  rotation  of  the  scene.  This  is  a  hololens  contribution  -  not  due  to  the 
scene  itself.  This  had  not  been  observed  before  and  required  some  investiga¬ 
tion.  It  was  observed  that  this  background  was  always  much  larger  when  all  25 
elements  were  exposed  simultaneously  with  the  same  input  scene  and  the  filters 
thus  made  addressed  in  parallel.  This  is  quite  likely  a  coherence  effect.  The 
25  separate  beams  from  the  filters  may  add  somewhat  coherently  at  the  detec¬ 
tor  thus  producing  a  larger  background  intensity  than  would  be  expected. 

It  was  known  that  the  correlation  spot  observed  on  the  television  monitor  was 
translation-insensitive  to  a  large  degree,  that  is,  the  input  scene  could  be 
linearly  translated  without  losing  the  correlation  spot.  The  bright  spot  would 
move  along  with  the  scene.  The  first  test  was  to  observe  whether  the  spot 
would  still  track  the  scene.  When  the  scene  was  translated,  an  interesting 
result  took  place,  which  perhaps  should  have  been  expected:  the  correlations 
separated.  The  bright  spot  contribution  due  to  the  scene  moved  as  always,  but 
another  portion  of  the  spot  did  not.  This  is  illustrated  in  Figure  35.  The 
scene  correlation  still  possessed  the  normal  rotational  variance;  but  the 
"noise"  correlation  did  not.  It  remained  fairly  constant  through  a  360  degree 
rotation  of  the  scene.  It  was  concluded  that  the  noise  correlation  was  due  to 
the  holoiens  itself.  Thus,  the  correlation  intensity  observed  on  the  televi¬ 
sion  monitor  consists  of  two  separable  parts  -  that  due  to  the  scene  itself  and 
that  due  to  the  hololens.  It  was  believed  that  any  imperfections  (dirt, 
smudges  or  artifacts  from  the  production)  would  show  up  in  the  correlation 
plane  as  a  background.  This  contribution  may  be  very  small  when  only  one  ele¬ 
ment  of  the  5X5  array  is  used  but  it  becomes  quite  large  when  aii  25  elements 
are  employed.  In  order  to  test  the  hololens  noise  contribution  idea,  it  was 
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Figure 


34 .  The  Ci  n-iation  intensity  from  all  25  elements  versus  the  angular 
rotati  'ii  of  the  input  scene.  The  scene  correlation  and  background 
correlation  due  to  the  hololens  are  superimposed. 
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necessary  to  try  to  duplicate  the  phenomenon  using  a  "noiseless"  lens.  This  was 
done  by  replacing  the  hololens  with  a  high  quality  glass  lens  of  approximately 
the  same  focal  length  and  adding  "noise"  to  the  clean  lens.  The  added  noise 
consisted  of  three  regular  array  mesh  screens  placed  directly  in  contact  with 
the  front  surface  of  the  lens  -  one  at  a  time.  Each  of  the  screens  contained  a 
different  fundamental  spatial  frequency.  It  was  observed  that  all  three 
screens  produced  background  correlations  very  similar  to  those  observed  using 
the  hololens.  The  highest  spatial  frequency  screen  seemed  to  produce  the 
largest,  most  localized  background.  Perhaps  this  indicates  that  the  noise  con¬ 
tained  in  the  hololens  was  also  of  a  high  spatial  frequency  content. 

An  interesting  pictorial  result  of  using  the  25  element  array  is 
given  in  Figure  36.  The  input  scene  was  a  military  tank  (side  view).  The  pho¬ 
tograph  was  made  by  placing  a  Polariod  (type  52)  film  a  few  inches  after  (for 
magnification)  the  next  image  plane  past  the  Fourier  transform  plane  of  the 
hololens  and  exposing  for  about  .2  seconds.  This  graphically  illustrates  the 
hololens's  ability  to  produce  25  images  from  a  single  input  scene. 

F.  Twenty-Five  Different  Filters 

The  obvious  next  step  was  to  attempt  to  store  25  different  filters 
and  address  all  of  them  simultaneously.  This  was  accomplished  after  experimen¬ 
tation  with  the  computer-controlled  rotation  stage  and  exposure  times.  The 
input  scene  in  this  instance  was  the  aerial  view  of  Huntsville.  Exposures  were 
made  for  every  14  degrees  of  rotation  of  the  scene.  The  K  ratio  was  again 
about  1  and  the  exposure  time  a  constant  12  seconds.  Only  one  element  of  the 
5X5  array  was  used  per  exposure.  This  was  controlled  by  the  mask/translation 
stage  combination  described  earlier.  The  first  experiment  was  to  access  all  25 
filters  simultaneously,  but  a  problem  arose.  Many  of  the  correlations  were 
excessively  weak  and  completely  masked  by  electrical  and  hololens  noise.  A 
translation  of  the  input  scene  might  have  alleviated  some  of  the  problem 
(hololens),  but  was  not  done.  It  was  decided,  rather,  to  access  the  filters 
one  at  a  time  (series  access)  and  try  to  obtain  better  correlations.  The 
results  are  given  in  Figure  37.  As  one  can  see,  the  peaks  are  clearly  visible. 
The  hololens  noise  is  minimal  when  the  filters  are  accessed  in  series.  It  was 
believed  at  this  stage  of  the  investigation  that  the  filters  might  be  too 
closely  spaced  (in  rotation  of  the  scene)  for  one  filter  to  be  distinguished 
from  another.  This  proved  not  to  be  the  case,  as  will  be  discussed  in  a  later 
section.  It  was  observed  that  the  correlation  curves  of  Figure  37  had  a 
smaller  FWHM  than  previously  measured.  This  is  probably  due  to  a  small  change 
in  scale  of  the  input  scene  due  to  refocusing  and  adjustment  to  the  zoom  lens 
of  the  input  camera.  The  correlation  as  a  function  of  scale  change  was 
discussed  in  Section  II  and  will  be  further  discussed  later  in  this  Section. 

The  important  aspect  is  that  25  filters  were  stored  on  the  glass  plate  in  an 
area  of  about  1.2  cm2  _  iess  than  the  area  of  a  thumbnail.  A  high  speed  method 
of  accessing  these  filters  in  series  would  make  even  this  early  result  useful 
for  a  variety  of  applications. 

The  goal,  of  course,  was  to  access  these  filters  in  parallel  -  not 
series.  This  was  attempted  using  the  same  scene  under  Identical  circumstances 
except  that  the  amount  of  scene  rotation  between  exposures  was  increased  from 
14  to  20  degrees  and  the  exposure  time  was  shortened  to  10  seconds  Instead  of 
12.  This  produced  the  parallel-addressed  correlations  of  Figure  38.  The  zero 
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correlation  intensity  axis  represents  the  background  due  to  the  hololens.  A 
possible  explanation  for  this  result  is  that  perhaps  the  shorter  exposure  time 
allowed  less  of  the  high  spatial  frequency  components  (associated  with  the 
hololens)  to  be  recorded  in  the  filter.  Even  though,  as  stated  earlier,  this 
background  is  fairly  constant,  even  small  variations  are  enough  to  destroy  an 
already  weak  correlation  signal. 

G.  Seventy-Three  Different  Filters 

The  success  of  the  25  filter  correlation  experiment  prompted  investi¬ 
gation  into  storing  even  more  information.  The  attempt  to  store  75  filters 
required  exposing  each  element  of  the  array  three  times.  Preliminary  investi¬ 
gation  of  this  procedure,  and  use  of  Equation  (84),  determined  that  the  expo¬ 
sure  times  should  be  3  seconds  for  the  first  exposure,  4  seconds  for  the  second 
and  5  seconds  for  the  third.  The  subject  was  a  scale  model  of  a  T62A  military 
tank.  The  average  power  in  each  of  the  25  elements  was  .12  W.  The  LCLV  was 
set  at  f  =  1400  Hz  @  6.5  volts,  as  usual.  The  optimum  K  value  was  experimen¬ 
tally  determined  to  be  about  1.0.  Exposures  were  made  for  every  4.6  degrees  of 
rotation  of  the  tank.  The  exposures,  rotation,  and  movement  of  the  photo¬ 
graphic  plate  were  computer-controlled  so  that  after  making  all  the  preliminary 
adjustments,  the  entire  75  exposures  took  only  about  15  minutes.  Results  of 
the  experiment  after  developing  and  placing  the  filter  back  in  the  correlator 
are  given  in  Figure  39.  All  25  elements  were  addressed  in  parallel.  Seventy- 
three  of  the  correlations  are  quite  distinct.  The  very  weak  correlations 
occurring  on  the  graph  are  those  corresponding  to  weak  elements  in  the  array. 

A  higher  quality  hololens  would  no  doubt  allow  all  75  filters  to  correlate.  It 
was  also  observed  that  the  background  did  not  overwhelm  the  correlations  as  had 
been  the  case  before  when  all  25  elements  were  addressed  in  parallel.  Perhaps 
the  lower  spatial  frequency  input  scene  caused  this  effect.  The  correlation  of 
more  than  70  matched  filters  on  a  single  plate  represents  a  vast  (fivefold) 
increase  over  any  similar  object  recognition  system  reported  to  date. 

H.  Sensitivity  to  Scale 

A  model  was  given  in  Section  II  for  the  effect  that  a  scale  change 
should  have  upon  the  correlation  intensity.  The  model  was  for  a  highly  spe¬ 
cialized  input  scene  (Ronchi  ruling).  In  order  to  test  this  model  using 
experimental  techniques,  a  Ronchi  ruling  type  input  scene  was  used  in  the 
correlator  and  matched  filters  made  as  before.  The  ruling  was  composed  of  flat 
black  strips  of  tape  on  a  white  paper  background.  The  black  strips  were 
one-quarter  inch  wide  with  a  one-quarter  inch  wide  separation  between  them. 

The  image  produced  by  the  LCLV  had  a  fundamental  spatial  frequency  of  .73  per 
mm  -that  is,  the  dark  areas  were  about  0.68  mm  wide  and  separated  by  the  same 
distance.  All  25  elements  were  exposed  simultaneously  as  in  Paragraph  E  of 
this  section.  The  matched  filters  made  looked  very  much  like  the  photograph  of 
Figure  29.  The  separation  between  the  zeroth  and  the  first  order  was  measured 
to  be  about  .013  cm.  An  exposure  time  of  12  seconds  with  K  *  1  was  found  to 
produce  good  correlations.  The  scale  (magnification)  of  the  input  scene  was 
then  changed  using  a  motorized  zoom  lens  made  by  C0HU  Corporation.  The  corre¬ 
lation  intensity  was  measured  using  a  Keithley  (Model  179)  digital  voltmeter 
connected  to  the  Colorado  Video  image  digitizer.  The  magnification  was  deter¬ 
mined  by  direct  measurement  of  the  input  scene  on  a  television  monitor.  Figure 
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40  gives  the  experimental  data  and  the  prediction  of  the  model  given  in  Section 
II.  The  agreement  is  better  than  expected  given  the  crudeness  of  the  model, 
which  neglects  all  but  the  fundamental  spatial  frequency.  Scale  changes  for  a 
more  "real  world"  type  scene  were  done  later  and  those  results  are  given  in 
Figure  41.  The  general  trands  are  similar.  Of  course,  the  widths  of  the 
curves  are  different,  the  Huntsville  scene  having  a  higher  spatial  frequency 
content  than  the  Ronchi  ruling.  The  downward  pointing  arrows  on  the  x  axis 
indicate  that  the  true  zero  is  somewhat  below  the  axis  in  Figure  34.  This  is 
due  to  the  hololens  noise  previously  discussed.  Notice  that  a  scale  change  of 
only  2  percent  is  necessary  to  decrease  the  correlation  by  half  for  the 
Huntsville  scene. 

I.  Sensitivity  to  Tilt 

The  previous  experiment  was  repeated  with  tilt  as  the  variable.  The 
same  Ronchi  ruling  input  was  used  as  before.  The  tilt  model  of  Section  II 
proved  again  to  be  reasonably  valid  for  this  specialized  input  scene.  Tilting 
the  input  scene  was  attempted  using  the  Scientific  Atlanta  stage  descibed  in 
Section  III,  but  a  problem  arose.  Tilting  the  scene  with  this  stage  is 
accomplished  via  a  pivot  below  the  stage.  Thus,  when  the  stage  is  tilted,  it 
is  also  translated  and  this  should  not  be  allowed  for  the  simple  model  to  be 
valid.  This  was  corrected  by  using  another  stage  which  had  its  pivot  in  the 
same  plane  as  the  stage.  Unfortunately,  this  stage  was  not  automated  so  data 
had  to  be  taken  manually  by  reading  the  voltage,  normally  going  to  the  chart 
recorder,  on  a  digital  voltmeter.  This  proved  to  be  adequate  even  though  con¬ 
tinuous  curves  were  not  possible.  Figure  42  gives  the  results  for  the  Ronchi 
ruling  input  scene;  again  the  fit  is  marginally  good.  Figure  43  gives  the 
results  for  the  Huntsville  scene  discussed  previously.  The  FWHM  for  the 
Huntsville  correlation  is  about  16  degrees  -  slightly  less  than  that  for  the 
Ronchi  ruling  scene,  as  should  be  expected. 

J.  Sensitivity  to  Rotation 

Throughout  this  investigation,  most  of  the  data  has  been  taken  with 
rotation  of  the  input  scene  being  the  variable.  The  simple  model  of  Section  II 
can  be  used  to  predict  the  change  in  correlation  intensity  with  rotation  of  the 
scene  when  the  scene  is  fairly  specialized  as  in  the  preceding  sections  of  this 
section.  Figure  44  compares  the  theory  to  experimental  results  for  a  Ronchi 
ruling  type  input  scene.  The  results  are  similar  to  those  obtained  in  the  pre¬ 
ceding  two  sections  except  for  the  widths  of  the  curves.  In  this  experiment  it 
was  possible  to  use  the  motorized  rotation  stage,  although  experimental  data  is 
plotted  as  points  to  be  consistent  with  the  previous  two  sections.  Figure  45 
gives  the  results  for  the  experiment  performed  using  the  Huntsville  scene.  The 
FWHM  is  observed  to  be  about  2.9  degrees  which  is  consistent  with  that  of 
Figure  32  -  a  correlation  produced  by  a  different  filter  of  the  same  scene.  It 
is  interesting  to  compare  results  from  this  section  to  those  of  the  previous. 
The  two  experiments  were  done  as  closely  alike  as  possible  in  K  ratio,  exposure 
times,  etc.  This  shows  that  the  correlation  is  much  less  sensitive  to  tilt 
angle  than  to  rotation.  This  is  an  important  result,  because  it  means  that 
fewer  matched  filters  are  needed  to  keep  an  objective  "in  sight"  when  only  the 
angle  of  approach  is  important.  On  the  other  hand,  it  also  means  the  angle  of 
approach  cannot  be  known  with  a  very  great  accuracy  in  actual  circumstances. 


Figure  41.  Correlation  versus  scale  change  (magnification)  for  a 
filter  of  the  aerial  map  of  Huntsville. 
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Figure  44.  Theoretical  prediction  and  experimental  results  —  correlation 
versus  rotation  for  a  Ronchi  ruling  input  scene. 


Correlation  versus  rotation  for  a  matched  filter  made 
of  the  aerial  Huntsville  ,  AL,  map. 


K.  Resolving  Closely  Spaced  Optical  Correlations 

The  ability  to  store  a  large  number  of  matched  filters  and  obtain 
correlations  of  them  presented  a  problem  which  had  not  been  encountered 
previously.  If  many  matched  filters  are  made  of  a  scene  only  slightly  rotated 
after  each  exposure  -  how  are  the  correlations  to  be  resolved  from  each  other? 
Often  the  orientation  of  the  scene  is  of  primary  interest.  Many  filters  must 
often  be  used  so  that  the  recognition  device  never  completely  "loses  sight”  of 
the  target  for  any  long  period  of  time.  The  experiment  consisted  of  doubly- 
exposing  a  photographic  plate  (overlapping  matched  filters)  to  an  input  scene 
at  two  slightly  different  rotations  positions.  Classic  lineshapes  were  used  to 
model  examples  of  experimental  results  and  some  conclusions  are  drawn  using  the 
Rayleigh  criterion  as  a  guideline. 

A  Gaussian  intensity  distribution  has  the  form  [58] 


II  =  Ic  e 


o  2/  2 
-2r  /w 


(98) 


where  r  is  the  radius  from  the  center  of  the  distribution  and  2w  is  the  full 
linewidth  at  I  =  l/e2  IQ.  If  another  distribution  of  the  same  form  is  located 
a  small  distance  away 
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then  the  sum  of  these  distributions  is 
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Pl  is  thus  the  ratio  of  the  peak  at  r  =  0  to  the  height  of  the  valley  between 
the  two  distributions.  Solving  for  ,  using  the  fact  that  the  second  term  is 
much  smaller  than  the  first,  gives 


6j  =  w  [2  Ln  (2p)]V2 


(103) 


Using  the  Rayleigh  criterion  as  a  rough  guideline,  pj  =  1.23,  [59]  and  substi¬ 
tuting  this  into  Equation  (6)  yields 
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(104) 
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The  same  sort  of  analysis  may  be  applied  to  a  Lorentzian  distribution  having 
the  form  [58]  (p.  154) 


I  = 


r2  +  (o/2)2 
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(105) 


where  A  is  a  constant  equal  to  I0(o/2)2,  r  is  the  radius  and  a  is  the  full 
linewidth  at  half  the  maximum  intensity.  This  produces  a  &2: 
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Then  for  P2  =  1.23  as  a  value  of  the  intensity  ratio 
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For  a  sinc^  distribution  [60] 
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where,  again,  r  is  the  radius  from  the  center  of  the  distribution  and 

a  =  x/r(o)  (109) 

where  r(o)  is  the  value  of  r  at  the  first  zero  of  intensity.  This  distribution 
yields 


a  <5 +  sin^T 

p3  =  - 1 i  (110) 

8  8in2  (063/2) 

which  must  be  solved  graphically  or  by  iteration  using  P3  *  1.23.  For  small 
values  of  63,  note  that  a 2  62  »  sin263 


Figure  46  shows  these  three  distributions  compared  with  experimental 
data  taken  as  before.  All  three  distributions  have  been  arbitrarily  fit  to 
experimental  data  at  r  *  0  and  r  »  2.95°  using  w  =  2.80°,  o  *  2.07°,  and  a  * 
.775  deg-1.  This  gives  6j  *  3.8°,  62  =  2.13°,  and  63  -  4.0°.  Note  that 
these  widths  should  depend  upon  the  spatial  frequency  of  the  input  scene  as 
well  as  other  factors  and,  thus,  serve  only  as  examples  here.  In  this  par¬ 
ticular  experiment  the  input  scene  was  a  low  resolution  black  and  white  aerial 
photograph  of  Huntsville,  AL.  The  photograph  was  placed  on  a  rotatable  stage 
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Figure  46.  Correlation  versus  rotation  for  Huntsville  aerial  map 
input  scene  compared  to  three  classical  lineshapes. 
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and  matched  filters  made  at  the  desired  rotation  angles.  Experimental  data  in 
Figure  46  is  the  correlation  intensity  for  a  filter  made  at  0  degree  rotation 
of  the  input  scene.  After  development  of  the  film  plate  (Kodak  649F)  it  was 
replaced  in  the  correlator  and  the  input  scene  rotated  slowly  as  the 
corresponding  correlation  intensity  was  measured. 

In  order  to  investigate  multiply-stored  filters,  a  single  plate  was 
exposed  twice  with  different  rotation  angles  of  the  input  scene.  After 
development,  the  filter  was  reinserted  into  the  correlator  and  data  taken  as 
before.  Figure  47  shows  the  actual  data  and  a  plot  of  Equation  (105)  using  the 
Rayleigh  criterion.  The  exposures  were  made  using  the  same  criteria  for  62  of 
2.13  degrees  of  scene  rotation.  The  agreement  is  reasonably  good.  Some 
variation  is  expected  due  to  the  difficulty  in  obtaining  equally  intense  corre¬ 
lations  for  the  two  filters.  This  depends  upon  the  exposure  times  used,  film 
response,  and  several  other  variables.  In  the  theoretical  analysis  the  calcu¬ 
lations  could  easily  have  been  done  using  distributions  having  different  peak 
heights  but  then  the  value  of  6  would  depend  upon  these  heights  which  cannot  be 
known  with  any  certainty  in  advance  of  doing  the  experiment.  In  the  calcula¬ 
tion  presented,  it  was  also  assumed  that  individual  correlation  intensities 
should  be  added  directly  rather  that  adding  the  corresponsing  fields  then 
squaring  for  the  intensity.  This  has  been  justified  experimentally  by  storing 
one  filter  made  at  0  degree  scene  rotation  and  another  filter  at  a  90  degree 
scene  rotation,  then  adding  the  resulting  correlation  intensities 
algebraically.  This  distribution  was  then  compared  with  data  similar  to  Figure 
47.  The  results  were  found  to  be  essentially  equal.  This  indicates  that  there 
is  iittie,  if  any,  phase  addition  contribution  to  the  sum  of  the  two  correla¬ 
tion  intensities  from  multiply  exposed  filters. 

In  this  section,  results  have  been  presented  from  several  experiments 
dealing  with  the  production  and  analysis  of  optical  matched  filters  made  using 
a  holographic  lens.  The  hololens  Itself  was  tested  to  determine  If  it  could 
reproduce  well-known  Fourier  transforms  and  the  efficiency  of  the  lens  was 
measured.  The  matched  filters  produced  by  the  hololens  were  then  compared  with 
those  made  using  a  conventional  glass  lens.  Finally,  the  hololens  was  used  in 
a  variety  of  experiments  that  until  now  could  not  have  been  done.  Large  num¬ 
bers  of  matched  filters  were  stored  on  photographic  plates  and  addressed  in 
parallel  and  series.  These  filters  produced  correlations  that  were  tested  for 
sensitivity  to  rotation,  scale,  and  tilt  of  the  input  scene.  Matched  filters 
were  made  of  a  simple  Ronchi  ruling  and  compared  to  the  theoretical  predictions 
of  Section  II.  Lastly,  an  investigation  was  made  Into  resolving  correlations 
produced  by  scenes  separated  only  a  few  degrees  in  rotaion.  It  was  found  that 
2.13  degrees  of  rotation  of  the  Huntsvllie  scene  was  necessary  to  satisfy  the 
Rayleigh  criterion  and  thus  produce  easily  separable  correlations. 
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DISCUSSION 


A.  Theory 

A  discussion  of  the  work  presented  here  should  begin  with  the  theoret¬ 
ical  anaiysis  used  to  modei  some  of  the  results.  The  idea  for  the  overlapping 
circles  came  from  a  paper  by  Casasent  [42].  In  that  paper  the  model  was  used 
to  predict  the  correlation  intensity  decrease  as  the  input  scene  was  rotated. 
Unfortunately,  the  model  was  fit  to  the  correlation  of  a  complex,  real  world 
scene  -  the  situation  where  the  modei  should  be  least  applicable.  Imagining 
that  the  Fourier  transform  of  a  simple  Ronchi  ruling  consists  of  circular  spots 
might  not  be  so  unreasonable  but  the  Fourier  transform  of  a  complicated  scene 
such  as  the  aerial  view  of  Huntsville  does  not  resemble  a  series  of  circular 
spots  by  any  means.  It  is  for  this  reason  that  none  of  the  correlations  of 
complicated  scenes  were  modeled  in  this  investigation.  A  model  which  has  the 
capability  of  handling  hundreds  of  spatial  frequencies  is  probably  achievable 
but  beyond  the  scope  of  this  work.  The  simple  model  of  Casasent  was  quite  ade¬ 
quate  in  predicting  general  trends,  however.  In  this  work,  the  rotational 
model  of  Casasent 's  was  expanded  so  as  to  predict  scale  and  tilt  variance  in 
correlation  and  gave  some  insight  into  the  results  obtained  when  complex  scenes 
were  used. 

There  are  several  reasons  for  performing  the  correlation  calculations 
in  the  Fourier  transform  plane  instead  of  the  image  plane.  The  correlation 
actually  occurred  in  the  Fourier  transform  plane,  that  is,  matched  filters  were 
made  of  the  Fourier  transform  of  the  scene,  not  the  scene  itself.  The  calcula¬ 
tions  were  also  much  easier  to  do  in  the  Fourier  transform  plane.  The  scale 
change  and  tilt  calculations  would  be  especially  cumbersome  to  do  in  the  image 
plane,  while  in  the  Fourier  transform  plane  these  effects  are  viewed  as  a 
change  in  location  of  the  spatial  frequency  components  characteristics  of  the 
image.  The  results  should  be  the  same  whether  the  calculations  are  done  for 
the  image  plane  or  the  Fourier  transform  plane.  This  was  found  to  be  true  for 
the  rotation  model  when  the  calculations  were  done  both  ways. 

B.  Holographic  Lens  and  LCLV 

This  work  represents  the  first  attempt  to  study  a  multiple  focus 
holographic  lens  In  a  real  time  optical  correlator.  The  results  of  tests  per¬ 
formed  on  the  hololens  showed  that  it  could  perform  almost  as  well  as  a  conven¬ 
tional  glass  lens  in  a  similar  situation.  The  hololens  faithfully  produced  the 
Fourier  transform  of  well-known  input  scenes  such  as  a  circular  aperture  and  a 
Ronchi  ruling.  It  was  encouraging  to  find  that  measurements  taken  from  the 
photograph  of  the  Fourier  transform  agreed  with  the  calculations  to  within  about 
10  percent.  This  10  percent  could  easily  have  come  from  errors  associated  with 
determining  the  magnification  of  the  input  scene  before  it  addressed  the  holo¬ 
lens  or  In  the  direct  measurement  of  the  separation  of  the  rings  in  the  Airy 
disk  (or  lobes  produced  by  the  Ronchi  ruling). 

The  modulation  transfer  function  (MTF)  of  the  hololens  was  not 
measured  directly.  In  order  to  do  this  type  of  measurement,  an  input  scene 
with  a  higher  resolution  than  that  of  the  television  monitor  and  LCLV  would  be 
needed.  In  routine  observations  of  the  image  produced  by  the  LCLV  it  seemed 


that  the  contrast  versus  spatial  frequency  (the  MTF)  [49]  (p.  118),  was  limited 
by  the  television  monitor,  not  the  LCLV.  Consider  a  typical  television  monitor 
measuring  about  10-inches  across  and  containing  512  lines  alternating  black  and 
white.  The  maximum  fundamental  spatial  frequency  produced  is  about  2  per  mm. 
When  focused  onto  a  1-inch  LCLV  this  produces  a  fundamental  spatial  frequency 
of  about  20  per  ram  which  is  well  within  the  capability  of  the  quoted  resolution 
of  the  LCLV.  The  contrast  at  this  spatial  frequency  is  not  known  precisely  for 
the  LCLV,  but  it  is  better  than  50  percent  [28].  Bleha  quotes  the  liquid 
crystal  light  valve  he  tested  (very  similar  to  the  one  used  in  this  work)  as 
having  60  lines/mm  resolution  at  50  percent  MTF  so  at  about  20  lines/mm  the  MTF 
is  prob-ably  much  higher  than  50  percent. 

Experiments  reported  in  Section  IV  of  the  present  investigation  using 
a  Ronchi  ruling  input  scene  did  demonstrate  some  of  the  capabilities  of  the 
LCLV.  The  LCLV  could  faithfully  produce  an  image  having  a  fundamental  spatial 
frequency  of  about  6.3  per  cm  (0.63  per  ram)  and  maintain  a  very  sharp  contrast. 
This  spatial  frequency  was  then  increased  by  zooming  the  camera  lens  away  from 
the  Ronchi  ruling  scene  until  the  fundamental  spatial  frequency  addressing  the 
LCLV  was  slightly  less  than  20  per  mm  (about  2  per  mm  on  the  TV  monitor).  A 
Ronchi  ruling  having  a  fundamental  spatial  frequency  of  2  per  mm  was  barely 
distinguishable  on  the  television  monitor  and  the  image  was  of  poor  contrast. 
The  image  produced  by  the  LCLV  was  also  barely  distinguishable.  This  investi¬ 
gation  indicated  that  the  MTF  of  the  LCLV  is  at  least  as  good  as  the  television 
monitor,  if  not  better. 

The  MTF  of  the  hololens  was  examined  in  a  similar  manner.  It  was 
found  that  by  observing  the  image  two  focal  lengths  away  from  the  lens  (an 
image  plane)  the  hololens  could  reproduce  anything  the  television  monitor/LCLV 
combination  was  capable  of  sending  it.  All  25  images  were  good  (see  Figure  36) 
with  the  overall  intensities  of  some  being  less  than  others  as  would  be 
expected  from  the  testing  reported  in  Section  IV.  While  the  MTF  of  the  optical 
system  was  not  measured  in  a  very  systematic  manner  it  can  be  said  that  the  MTF 
of  the  hololens  and  LCLV  was  probably  at  least  as  good  as  the  television  camera 
and  monitor  used  as  the  input. 

The  efficiency  and  intensity  distribution  produced  by  the  hololens 
proved  to  be  adequate  but  considerable  improvement  should  be  expected  in  the 
next  generation  of  hololens.  Discussions  with  the  manufacturer  have  led  to 
placing  an  order  for  a  hololens  which  produces  a  7X7  array  with  an  overall 
efficiency  better  than  15  percent  and  a  more  uniform  intensity  distribution 
[61].  This  should  lead  to  about  150  stored  matched  filters.  This  is  enough 
matched  filters  to  ensure  no  more  than  a  50  percent  decrease  in  correlation 
intensity  regardless  of  the  rotational  position  of  the  scene  with  the  present 
optical  correlator  and  a  scene  similar  to  the  aerial  view  of  Huntsville.  The 
FWHM  of  the  correlation  intensity  versus  rotation  curve  for  the  Huntsville 
scene  was  about  2.9  degrees  (Figure  45)  so  that  126  filters  would  be  needed  for 
no  more  than  50  percent  decrease  in  correlation  intensity.  For  a  tilt  of  _+  60 
degrees,  about  57  filters  would  be  needed  and  about  25  filters  would  be 
required  for  a  100  percent  change  in  scale.  This,  unfortunately,  does  not  tell 
the  whole  story.  Suppose  that  at  each  rotational  position  it  is  desired  that 
the  correlation  intensity  be  tilt-invariant.  A  total  of  more  than  7000  filters 


would  be  needed.  When  scale  change  is  Included,  the  problem  is  magnified  even 
more.  For  very  specialized  applications,  the  147  filters  may  be  enough, 
however.  In  fact,  for  a  stationary  aerial  observation  the  75  filters  investi¬ 
gated  here  would  produce  about  a  10  to  15  percent  correlation,  which,  for  some 
applications  might  be  sufficient. 

The  background  correlation  signal  observed  using  the  hololens  did  not 
prove  to  be  a  major  problem  in  this  investigation.  Fortunately,  this 
background  was  fairly  constant  throughout  a  full  360  degree  rotation  and  the 
scene  correlations  were  usually  distinguishable  above  the  background.  However, 
if  this  were  not  present,  perhaps  more  than  three  filters  could  be  stored  at 
each  location.  This,  of  course,  would  increase  the  overall  recognition  capabil¬ 
ities  of  the  system.  More  than  three  filters  at  each  location  were  attempted, 
but  fewer  than  half  the  correlations  were  of  significant  intensity. 

C.  Matched  Filters 

The  experimental  results  obtained  in  this  research  covered  several 
aspects  of  matched  filtering.  The  results,  however,  are  much  more  graphic  than 
perhaps  has  been  indicated  in  this  report.  In  actual  practice  the  observation 
of  the  correlation  spot  on  the  television  monitor  is  quite  dramatic.  The 
correlation  intensity's  sensitivity  to  rotation,  tilt  and  scale  is  extremely 
pronounced.  For  this  reason,  many  of  the  experimental  results  were  videotaped. 

The  making  of  the  initial  matched  filters  in  this  investigation 
required  a  great  deal  of  trial  and  error  with  exposure  times  and  developing 
techniques.  Fortunately,  many  of  the  procedures  were  known  from  earlier  work 
done  by  laboratory  personnel  with  conventional  photography. 

Early  in  the  research  it  was  obvious  that  the  correlations  obtained 
using  the  present  arrangement  were  different  in  some  respects  from  those 
obtained  by  the  Harris  group  [24].  A  major  difference  was  the  correlation 
widths.  The  Harris  group,  using  a  similar  input  scene,  obtained  widths  much 
narrower  than  those  presented  here  (fractions  of  a  degree)  for  rotating  the 
input  scene.  The  major  difference,  it  turns  out,  was  the  fact  that  that  group 
used  high  resolution  35  mm  slides  as  the  input  scene.  This  enabled  them  to 
record  much  higher  spatial  frequencies  than  the  commercial  television  monitor 
used  in  this  work.  Section  II  discussed  how  this  should  affect  the  correla¬ 
tions  obtained;  the  higher  the  spatial  frequency,  the  narrower  the  width.  The 
use  of  a  television  monitor,  however,  has  its  advantages.  The  input  scene  can 
be  instantly  changed.  The  input  scene  can  be  a  moving  object  and,  of  course, 
using  the  television  monitor  is  much  more  convenient  than  having  to  produce  a 
35  mm  slide.  The  television  monitor  in  conjuction  with  the  LCLV  is  what  makes 
this  a  "real  time"  recognition  system. 

D.  Optical  Memory 

In  the  course  of  this  investigation,  it  was  possible  to  significantly 
increase  the  optical  memory  of  the  correlator  over  anything  attempted  in  the 
past.  Until  the  advent  of  the  hololens  the  only  way  to  achieve  more  than  one 
memory  was  through  multiply-exposing  the  photographic  plate.  This  was  limited 
to  about  eight  images  being  stored,  which  in  itself  is  quite  an  accomplishment. 
These  filters,  however,  could  only  be  correlated  with  an  input  scene  one  at  a 


time;  and  it  was  not  possible  to  determine  which  filter  was  being  correlated. 
The  hoioiens  arrangement,  as  described  in  this  work,  allowed  25  filters  to  be 
addressed  simultaneously  and  allowed  the  observer  to  determine,  visually  or 
electronically,  which  filter  was  producing  the  correlation.  For  the  first 
time  a  target  could  be  identified  by  observing  the  location  of  the  bright 
correlation  spot  alone.  This  was  accomplished  by  a  slight  defocusing  of  the 
final  lens  in  front  of  the  CCD  detector.  This  produced  a  spatial  distribution 
of  the  correlations  as  they  occurred  (see  Figure  48).  In  many  applications, 
however,  it  is  not  necessary  to  know  which  target  is  being  observed  -  only 
that  one  of  the  preselected  targets  is  being  observed  by  the  correlator. 
Adjusting  the  position  of  the  final  iens  in  front  of  the  CCD  detector  can 
bring  ail  25  possible  correlation  spots  to  the  same  location;  thus,  only  one 
detector  need  be  used  in  this  application. 

Expanding  the  memory  further  by  using  multiple  exposure  techniques 
in  conjunction  with  the  hoioiens  allows  the  total  number  of  recognizable 
scenes  to  be  increased  at  the  expense  of  not  being  able  to  tell  exactly  which 
scene  produced  the  correlation.  In  this  Investigation,  75  scenes  were 
attempted  by  storing  three  at  each  of  the  25  locations  produced  by  the 
hoioiens;  73  of  the  filters  produced  distinguishable  correlations.  The 
correlations  were  detected  with  the  CCD  iens  adjusted  so  that  all  correlation 
spots  would  occur  at  the  same  coordinates  on  the  TV  monitor.  If  this  lens  had 
been  repositioned,  it  would  have  been  possible  to  obtain  results  similar  to 
those  when  only  25  filters  were  stored  -  however,  the  exact  scene  being 
observed  would  not  have  been  identifiable.  It  would  have  been  possible, 
however,  to  narrow  the  possibilities  to  three  scenes.  The  7X7  hoioiens,  which 
is  being  produced,  should  allow  49  targets  to  be  distinguished  and  up  to  147 
or  more  total  scenes  to  be  stored,  without  absolute  dlstlnguishability. 

E.  Rotation,  Scale,  and  Tilt 

The  experiments  using  a  Ronchi  ruling  input  scene  were  suggested  as 
a  test  of  the  simple  theory  given  in  Section  II  [62],  Results  were  actually 
better  than  expected  given  the  simple  nature  of  the  model.  A  Ronchi  ruling 
had  been  used  very  early  in  the  investigation  to  test  the  hoioiens.  It  is 
realized  that  the  three  variables  (rotation,  scale,  and  tilt)  do  not  form  a 
complete  picture.  In  the  rotation  experiment,  the  other  two  variables  were 
kept  constant  -  similarly  for  the  scale  and  tilt  experiments.  Using  the 
experimental  results  obtained,  it  may  be  possible  to  Infer  the  effects  of  a 
simultaneous  rotation  and  scale  change  or  other  combinations  of  the  three 
variables.  It  was  shown  In  Section  IV  (Paragraph  K)  that  the  correlation 
intensities  could  be  added  directly.  Using  this  as  a  guideline,  it  should  be 
possible  to  predict  the  effect  of  rotation,  scale,  and  tilt  performed 
simultaneously.  As  an  example,  consider  the  Huntsville  scene  correlation 
results  of  Section  IV.  A  rotation  of  1  degree  might  produce  a  20  percent  loss 
in  correlation  Intensity;  a  scale  change  of  1  percent  might  cause  another  30 
percent  loss;  while  a  tilt  of  4  degrees  might  cause  a  15  percent  loss. 
Combining  ail  three  losses  one  would  expect  to  achieve  a  48  percent  coi rela¬ 
tion  on  the  scene.  This  result  is  somewhat  conjecture  but  seems  reasonable. 
Experiments  of  this  type  were  not  performed  but  should  be  simple  enough  to  do 
In  the  future. 


F.  Resolution  Between  Correlations 

Later  in  the  investigation  when  multiple  exposure  and  arrays  of 
filters  were  being  made,  it  was  observed  that  a  limit  existed  for  being  able 
to  separate  the  correlation  profiles  of  scenes  separated  by  only  a  few  degrees 
of  rotation.  An  experiment  was,  thus,  devised  and  performed  to  determine  what 
this  limit  was  for  a  specific  scene.  Results  of  this  experiment  were  given  in 
Section  IV  and  in  a  paper  [57].  The  experiment  was  done  using  rotation  as  the 
only  variable,  but  the  analysis  technique  should  carry  over  into  scale  and 
tilt  experiments.  The  idea  to  use  the  Rayleigh  criterion  seemed  to  be  quite 
natural  after  observing  the  correlation  curves  produced  as  in  Figure  47.  It 
was  necessary,  in  order  to  model  the  results,  to  assume  some  type  distribution 
of  intensity  versus  rotation.  Three  well-known  distributions  were  fit  to  the 
experimental  data.  The  best  fitting  distribution  (Lorentzian)  was  then  used 
to  predict  the  sum  of  two  correlation  curves.  It  should  be  mentioned, 
however,  that  the  shape  may  not  be  the  best  fit  for  other  input  scenes.  This 
experiment  was  done  mostly  to  point  out  the  need  for  resolution  criteria;  none 
exists  at  this  time.  The  Rayleigh  criterion  is  well  known  and  accepted  and 
should  serve  as  an  adequate  beginning  for  resolving  matched  filters. 
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VI 


CONCLUSIONS  AND  SUGGESTIONS  FOR  FUTURE  INVESTIGATIONS 


This  investigation  represents  the  initial  research  into  large  opticai 
memories  consisting  of  arrays  of  photographically  stored  matched  filters.  A 
number  of  specific  conclusions  can  be  drawn  from  the  work  given  in  this 
report. 


A  prototype  holographic  eiement  which  produced  a  5X5  array  of  Fourier 
transforms  of  a  scene  was  tested  and  found  to  be  acceptable  for  simple  optical 
correlation  applications.  Limits  of  the  holographic  element  were  found, 
however,  in  that  the  eiement  itseif  contained  imperfections,  which,  at  times, 
made  the  correlation  signal  difficult  or  impossible  to  detect. 

A  definite  source  of  the  background  noise  observed  was  not  discovered 
but  several  possibilities  exist.  The  general  cleanliness  of  the  optical 
arrangement  used  to  produce  the  hololens  is  certainly  important.  Any 
fingerprints,  dust  or  contaminated  chemicals  used  could  leave  artifacts  on  the 
hololens  which  wouid  appear  in  the  correlation  signal.  The  neutral  density 
filters  used  in  the  aperture  of  Figure  21  to  equalize  the  intensities  of  the 
individual  elements  could  possibly  add  an  unwanted  background  if  they  con¬ 
tained  scratches,  fingerprints,  etc. 

Cross  taik  between  adjacent  filters  was  not  observed.  Initially  it  was 
feared  that  the  higher  order  diffraction  elements,  which  were  charcterlstic 
of  the  diffraction  screen  used  to  make  the  holoiens,  might  overlap,  thus, 
allowing  adjacent  filters  to  interfere.  The  hololens  used  in  this  investiga¬ 
tion  did  contain  higher  order  diffraction  elements  but  they  were  lower  in 
intensity,  by  at  least  an  order  of  magnitude,  than  the  25  primary  elements. 

The  normal  exposure  times  used  did  not  allow  these  higher  orders  to  even  be 
detected  on  the  photographic  piate.  A  new  holographic  element,  also  produced 
by  H.  K.  Liu,  will  be  tested  in  the  near  future.  This  new  eiement  produces  a 
7X7  array  of  Fourier  transform  elements  which  will  allow  up  to  49  scenes  to  be 
stored  with  discrimination  between  them  and  more  than  147  scenes  may  be 
recognized  without  absolute  discrimination.  If  the  background  noise  is  lower 
on  this  new  holographic  element,  it  may  be  possible  to  expose  each  element 
more  than  three  times,  thus,  allowing  an  even  greater  number  of  scenes  to  be 
recognized. 

The  use  of  the  LCLV  allowed  the  filters  to  be  addressed  in  real  time  and 
the  correlations  to  be  measured  directly  while  varying  the  Input  scene. 
Experiments  of  this  type  Included  changing  the  angular  and  linear  location  of 
the  input  scene  as  well  as  the  magnification  while  observing  the  change  in 
intensity  of  the  correlation.  It  was  found  that  the  correlation  Intensity  was 
much  more  sensitive  to  rotational  motion  of  the  input  scene  than  to  tilt.  The 
simple  models  developed  in  Section  II  proved  adequate  in  predicting  the  trends 
observed  experimentally. 

The  LCLV  proved  to  be  adequate  in  producing  coherent  Images  from  an 
incoherent  input  scene,  at  least  to  the  limitation  of  the  television  screen 
Input.  A  higher  resolution  Input  scene  might  show  the  limitations  of  the 
LCLV.  The  television  screen,  however,  probably  represents  the  type  of  input 
one  might  expect  to  be  used  in  an  actual  application  such  as  automatic  target 


acquisition.  Higher  resoiution  television  screens  are  avaiiabie  and  might 
approach  the  LCLV's  iimits,  but  none  were  used  in  this  investigation.  Research 
is  underway  at  this  time  to  improve  the  resoiution  as  well  as  the  reaction 
time  of  the  LCLV.  If  the  potential  users  (and  current  funders)  agree  that  the 
LCLV  is  the  best  approach,  then  it  is  very  likely  that  advances  will  be  made 
rapidly.  At  this  time,  the  LCLV  is  in  competition  with  other  devices  and 
techniques,  such  as  erasable  photoplastics  and  white  light  processing,  where  a 
coherent  monochromatic  source  is  not  necessary  [63]. 

The  resoiution  between  correlations  produced  by  very  similar  scenes  was 
investigated  out  of  necessity,  giving  some  insights  that  might  have  otherwise 
gone  unnoticed.  Matched  filters  of  very  similar  scenes  were  recorded  and  the 
correlations  observed  took  the  shape  of  the  Rayleigh  double  star  curve.  This 
observation,  thus,  led  to  the  development  of  a  simple  model  which  fit  experi¬ 
mental  data  reasonably  well.  This  model  can  probably  be  modified  to  predict 
the  resoiution  needed  as  a  function  of  the  spatial  frequency  of  the  input 
scene,  it  seems  that  this  would  be  quite  useful  in  the  future. 

The  hololens  approach  combined  with  the  LCLV  has  proven  itself  to  be 
quite  worthy  of  consideration,  even  in  its  present  state,  for  use  by  various 
agencies  of  the  government  for  target  acquisition,  discrimination,  and 
tracking.  The  present  investigation  has  shown  that  more  than  70  matched 
filters  can  be  stored  in  a  very  small  area  and  addressed  in  parallel.  This 
investigation  has  also  shown  that  by  storing  only  one  filter  per  array  location 
paraiieL  discrimination  can  be  achieved.  This  was  accomplished  by  a  slight 
defocuslng  of  the  lens  just  in  front  of  the  small  CCD  camera.  A  visual 
discription  is  given  by  the  photograph  of  Figure  48.  If  a  detector  were 
placed  at  each  location  of  a  bright  spot,  automatic  discrimination  could  be 
achieved  without  the  presence  of  an  observer.  If  the  array  of  Figure  48  were 
focused  to  a  point  by  the  CCD  camera  lens,  then  only  one  detector  would  be 
necessary;  but  the  discrimination  would  be  lost. 

Other  investigators,  notably  at  Grumman  [64],  are  at  this  time,  experi¬ 
menting  with  a  hololens  of  their  own  design.  Thus  far,  they  have  used  a  2X2 
array  with  multiple  exposures  to  achieve  better  than  16  matched  filter 
correlations.  The  5X5  array  Investigated  here  represents  a  better  than 
fourfold  increase  over  this  number. 

The  field  of  optical  computing  has  received  a  much  needed  boost  in  the 
past  year,  largely  due  to  the  promising  future  of  holographic  elements  and 
coherent  light  converters  such  as  those  Investigated  in  this  research. 

Optical  recognition,  as  presented  in  this  report  is  but  one  area  in  the  larger 
field  of  optical  computing;  but  It  is  the  first  area  that  has  been  advanced  to 
the  point  of  implementation  in  government  and  industrial  projects. 

In  summation,  the  specific  accomplishments  of  this  investigation  were: 

1.  The  first  testing  of  a  holographic  element  that  could  produce  a 
5X5  array  of  Fourier  transformed  Images. 

2.  The  first  real  time  addressing  of  73  matched  filters. 
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3.  The  first  true  discrimination  in  real  time  by  an  optical  corre¬ 
lator  -  an  individual  scene  was  identified  by  the  spatial  location  of  the 
correlation  signal. 

4.  The  first  investigation  of  the  effects  of  rotation,  scale,  and 
tilt  on  the  correlation  of  a  scene  in  real  time  using  a  liquid  crystal  light 
valve. 

5.  The  first  investigation  of  how  two  correlation  signals  may  be 
distinguished  from  each  other  when  their  respective  scenes  are  very  similar. 

6.  The  first  investigation  of  using  multiple  exposure  techniques  to 
produce  matched  filters  of  a  scene  Fourier  transformed  by  a  holographic 
element. 

7.  The  first  investigation  of  noise  associated  with  a  Fourier 
transforming  holographic  element. 

A  few,  somewhat  surprising,  results  were  obtained  in  this  investigation. 
The  separation  of  the  correlation  signal  from  the  hololens  noise  was  observed 
when  the  input  scene  was  translated  slightly.  The  correlation  signal 
displayed  the  expected  translation  invariance  while  the  noise  signal  did  not. 
This  background  was  found  to  be  quite  constant  throughout  a  360  degree  rota¬ 
tion  of  the  input  scene. 

Defocusing  the  camera  used  to  detect  the  correlations  also  provided  an 
unexpected  result.  This  allowed  a  spatial  distribution  of  the  correlation 
signals  so  that  the  correlation  produced  by  an  individual  filter  could  be 
distinguished,  thus,  making  absolute  discrimination  possible. 

Some  uncertainties  involved  in  this  investigation  might  also  be  worthy 
of  more  consideration: 

1.  The  Videotek  image  digitizer  used  was  not  calibrated  for  abso¬ 
lute  measure  of  the  correlation  intensities;  only  relative  measurements  were 
possible.  A  systematic  method  for  measuring  these  intensities  should  provide 
a  more  quantitative  measure  of  the  recognition  abilities  of  an  optical 
correlator. 

2.  The  display  produced  by  the  LCLV  was  not  consistent  on  a  day-to- 
day  basis  and  often  had  to  be  adjusted  before  an  experiment  could  be  done.  A 
method  of  monitoring  the  image  produced  by  the  LCLV  could  possibly  be  employed 
as  the  exposures  are  made.  This  could  be  done  by  using  a  beamsplitter  to 
remove  a  percentage  of  the  image  intensity  which  could  be  monitored  visually 
as  the  experiment  progresses. 

3.  The  output  power  from  the  HeNe  laser  used  to  produce  the 
coherent  image  was  periodically  checked  but  not  monitored  continuously.  A 
slight  decrease  in  power  could  give  erroneous  decrease  in  the  correlation 
intensities  observed.  A  beamsplitter  and  detector  placed  in  the  output  beam 
of  the  laser  would  alleviate  this  problem.  The  detector  output  voltage  could 
be  continuously  displayed  by  a  chart  recorder  as  the  experiment  progresses. 
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Suggestions  for  further  work  on  the  topics  discussed  in  this  report 
might  include: 


1.  Repeating  the  experiments  dealing  with  rotation,  scale  and  tiit 
using  a  higher  resolution  input  scene  which  would  take  full  advantage  of  the 
LCLV's  quoted  resolution  of  60  lines/mm.  A  significant  narrowing  of  curves 
such  as  Figure  45  might  be  expected;  thus,  making  for  a  more  accurate  measure¬ 
ment  of  the  position  of  the  scene. 


2.  The  correlation  intensity  separation  from  the  holoiens  noise  as 
given  in  Figure  35  should  possibly  be  Investigated  more  closely  as  a  method 
for  reducing  the  background  signal.  This  background  signal  might  possibly  be 
used  advantageously  in  that  it  provides  a  general  location  where  detectors 
should  be  placed  for  the  detection  of  the  expected  correlation  signal. 


3.  A  5X5  array  of  individual  detectors  placed  in  the  correlation 
piane  to  detect,  independently,  the  intensities  displayed  in  Figure  48  would 
provide  the  opportunity  to  measure  directly  the  detection  crosstalk  that  might 
be  experienced  in  an  actual  application  of  this  technology.  Excessive 
crosstalk  could  possibiy  produce  recognition  errors. 


4.  Matched  filters  using  dichromated  gelatin  photographic 
techniques  should  be  investigated  for  use  with  the  hoiolens.  The  efficiency 
of  dichromated  gelatin  is  much  greater  than  the  649F  plates  used  in  this 
investigation.  More  efficient  matched  filters  would  produce  larger,  more 
easily  detectable  correlation  signals. 


5.  A  method  for  electronically  subtracting  the  holoiens  background 
signal  should  be  investigated.  A  device  such  as  a  Quantex  image  digitizer 
might  be  employed  for  this  purpose  as  it  has  the  capability  of  subtracting  a 
stored  scene  from  a  displayed  scene.  An  electronic  subtraction  of  background 
might  allow  a  larger  number  of  multiple  exposures  to  be  made  than  was  done  in 
this  investigation. 
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The  Fourier  Transform  of  a  Ronchi  Ruliing 

In  Section  II  a  simpie  correlation  theory  was  developed  which  was  based 
upon  the  overlap  of  two  circles.  Later,  in  Section  IV,  this  simple  model  was 
used  to  predict  the  correlation  versus  several  variables  for  a  Ronchi  ruling 
input  scene.  The  rationale  was  that  the  Fourier  transform  of  a  Ronchi  ruling 
produced  a  roughly  circular  series  of  exposures  on  a  photographic  plate.  In 
this  appendix,  the  Fourier  transform  of  a  finite  Ronchi  ruling  is  calculated 
exactly.  Many  standard  textbooks  give  solutions  to  this  problem  for  a  ruling 
extending  an  infinite  distance  but  the  solution  to  a  finite  ruling  is  worth 
examining  because  of  the  application  to  some  of  the  experiments  presented  in 
Section  IV. 

Figure  A-l  gives  the  transmission  versus  £  for  a  Ronchi  ruling.  For  a 
ruLIng  extending  over  an  infinite  distance,  this  transmission  may  be  written: 


t<5)  -  £  rect 

n=-°° 

where,  in  general,  the  rect  function  is  written: 


rect  (Z) 


0  otherwise 


(112) 


A  finite  ruling  is  best  represented  as  a  convolution  of  a  rect  function  with  a 
comb  function  bounded  by  a  rect  function  window  -  illustrated  by  the  dashed 
Line  Ln  Figure  A-2. 


t(«  -  [rect  (I)  *  comb  (t)j  .  rect(£)  (113) 

where,  In  general,  the  comb  function  is  written: 


comb  (z)  =  6(z-n)  .  (114) 

n=-® 

It  is  desired  to  Illuminate  this  Ronchi  ruling  with  plane  wave  unit  amplitude 
light,  then  Fourier-transform  the  transmitted  field  us.  .g  a  lens.  If  the 
ruling  is  placed  in  the  front  focal  plane  of  the  lens  then  the  field  in  the 
back  focal  plane  will  be  proportional  to  the  Fourier  transform  of  Equation 
(113): 
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\  n=-°° 
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This  is  piotted  in  Figure  A-2.  No  "adjustable  parameters"  were  used  but 
rather  experimental  values  of  the  ruling  were  substituted  for  f,  a,  and  N. 
Results  are  quite  good  for  the  separation  between  peaks.  In  Section  IV,  the 
separation  was  found  to  be  about  1.3  X  10”2  cm  between  the  first  order  and  the 
zeroth.  This  agrees  very  well  with  the  prediction  of  Figure  42.  Only  the  odd 
orders  are  present  -  as  should  be  expected.  A  square  wave  may  be  written  as  a 
sum  of  sines  or  cosines  [65]  (p.  244): 

f(S)  -  k  6  +  -  cos  QZ  -  ~  cos  3Q£  +  cos  5QF,)  ,  (133) 

2  \  x  3tt  5h  / 

or 

f(0  -  ~  (l  +  -sln  QZ,  +  sin  3Qf,  +^sin  5QfJ  .  (134) 

Z  \  71  Jit  5tt  / 

In  either  representation,  only  the  odd  orders  of  Q  are  present. 

The  width  of  the  first  order  is  also  important.  In  the  model  presented 
in  Section  II  and  used  in  Section  IV  the  vaiue  of  R,  the  radius  of  the  exposed 
first  order  was  taken  directly  from  photographs.  In  the  analysis  presented  in 
this  appendix,  a  computed  value  can  be  taken  directly  from  Figure  A-2.  The 
FWHM  of  the  first  order  is  about  7.4  X  10~6  M.  Therefore  R  is  equal  to  3.7  X 
10~4  cm.  This  represents  an  error  of  about  14  percent  when  compared  to 
experimental  values  taken  in  Section  IV. 

It  was  possible  to  experimentally  verify  the  calculated  Fourier  trans¬ 
form  given  In  Figure  A-2.  The  Ronchi  ruling  described  in  Section  IV  was  again 
used  as  the  input  to  the  optical  system  of  Figure  24,  with  the  reference  beam 
blocked.  A  small  telescope  was  placed  at  the  location  usually  reserved  for 
the  film  holder.  The  overall  intensity  of  the  laser  ILght  was  decreased, 
using  neutral  density  filters,  to  the  point  where  the  optical  Fourier  trans¬ 
form  could  be  viewed  by  the  naked  eye.  The  resulting  Image  was  then  projected 
Into  a  small  television  camera  and  displayed  on  a  television  set.  The  image 
was  quite  similar  to  the  photograph  of  Figure  29.  In  order  to  obtain  an 
intensity  profile  of  the  image,  the  signal  from  the  television  camera  was 
input  to  a  Tektronix  oscilloscope  which  could  display  scans  of  single  televi¬ 
sion  lines.  The  oscilloscope  trace  was  then  photographed.  This  is  shown  in 
Figure  A-3.  The  zeroth  order  and  successive  odd  orders  are  present.  This 
compares  welt  with  the  calculation  done  earlier  in  this  appendix  and  plotted  in 
Figure  A-2. 
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The  previous  analysis  seems  to  somewhat  justify  the  simple  overlapping 
circles  model  presented  in  Section  II.  The  higher  orders  present  in  the  fre¬ 
quency  spectrum  have  been  ignored.  The  intensity  of  the  third  order  Is  quite 
weak  compared  to  the  first;  therefore,  the  neglection  of  higher  orders  is  pro 
bably  justifiable. 

The  analysis  presented  in  this  appendix  was  satisfying  in  that  the 
separation  between  peaks  as  well  as  the  width  of  the  first  order  peaks  was 
calculated  and  proved  to  be  comparable  to  experimental  values. 


Film  Characteristics  and  Developing  Techniques 


Almost  all  the  experimental  work  done  in  the  investigation  involved  pho¬ 
tographic  plates.  It  is  for  this  reason  that  this  appendix  is  included.  Some 
of  the  techniques  given  here  are  well  established;  others  were  developed  and 
tested  as  the  need  arose. 

A.  Kodak  694F  Film  [66] 

Kodak  649F  spectroscopic  plates  used  in  this  investigation  are  some 
of  the  highest  resolution  plates  available  -  3000  lines/mm.  It  is  a  fairly 
thick  (17  pm)  emulsion  on  a  .040-inch  thick  glass  plate.  The  granularity  is 
given  as  less  than  5  where  30  is  considered  coarse.  It  is  the  lowest  in  granu¬ 
larity  of  any  Kodak  film.  It  is  used  for  a  wide  variety  of  spectroscopic 
applications.  The  spectral  response  is  fairly  flat  from  425  nm  to  675  nm. 

The  MTF  of  649F  film  is  greater  than  .70  as  far  out  as  2000/am  spatial  fre¬ 
quency  [49]  (p.  572).  The  sensitivity  is  given  as  8X10~5  j/cra 2  in  the 
spectral  response  region. 

B.  Developing  [65] 

The  developer  used  throughout  this  investigation  was  Kodak  D-19. 

The  receipt  for  this  developer  is: 


Water  (50°C) 

500  cc 

Kodak  Elon 

2  g 

Sodium  Sulfide 

90  g 

Hydroquinone 

8  g 

Sodium  Carbonate 

52.5  g 

Potassium  Bromide 

5  g 

Cold  water  to  make 

1  liter 

The  principle  ingredient  in  Kodak  Elon  developing  agent  is  p-methylaminophenol 
sulfate.  The  proportions  of  Elon  and  hydroquinone  in  a  developer  will  vary 
according  to  the  film  being  developed.  Elon  alone  produces  a  low  contrast 
image  with  some  fog.  Hydroquinone  alone  produces  a  high  contrast  image  with  a 
high  fog  level.  The  stop  bath  used  was  Kodak  indicator  stop  bath.  The 
receipt  for  this  is: 

Water  1  liter 

Acetic  acid  (29%)  48  ml 

The  function  of  the  acidic  stop  bath  is  to  stop  the  developing  procedure 
caused  by  the  alkaline  components  of  the  developer.  The  fixer  used  was  Hunt's 
Aerofix.  The  receipt  for  this  is: 


Water  (50  °C)  600  ml 

Sodium  thiosulfate  360  g 

Ammonium  chloride  50  g 

Sodium  sulfide  15  g 

Acetic  acid  (28%)  48  ml 

Boric  acid  crystals  7.5  g 

Potassium  alum  15  g 

Cold  water  to  make  1  liter 


The  developing  procedure  was  as  follows:  All  chemicals  were  kept  at  20  °C  by 
using  a  temperature-controlled  water  baLh.  The  first  four  steps  were  carried 
out  In  complete  darkness. 


1.  Immersion  with  agitation  in  D-19  developer  for  6  minutes. 

2.  Immersion  with  agitation  in  stop  bath  for  3  minutes. 

3.  Immersion  with  agitation  in  fixer  for  1  minute. 

4.  Flowing  water  rinse  for  at  least  10  minutes. 

5.  Final  rinse  in  absolute  ethanol  for  at  least  5  minutes. 

6.  Air  dry  for  at  least  20  minutes. 

The  previous  procedure  seemed  to  give  consistently  good  matched  filters. 
Slight  variations  of  this  procedure  are  possible  without  causing  problems. 


B-4 


DISTRIBUTION 


No.  of 
Copies 


Commander 

US  Army  Research  Office 
ATTN:  AMXRO-PH,  Dr.  R.  Lontz 
P.  0.  Box  12211 

Research  Triangle  Park,  NC  27709 

US  Army  Research  and  Standardization  Group  (Europe) 

ATTN:  AMXSN-E-RX,  LTC  D.  R.  Reinhard 
Box  65 

FPO  New  York  09510 
Commander 

US  Army  Materiel  Development  and  Readiness  Command 
ATTN:  Dr.  James  Bender 

Dr.  Gordon  Bushey 
5001  Eisenhower  Avenue 
Alexandria,  VA  22333 

Headquarters 
Department  of  the  Army 

Office  of  the  DCS  for  Research,  Devleopment  &  Acquisition 

ATTN:  DAMA-ARZ 

Room  3A474 ,  The  Pentagon 

Washington,  DC  20301 

OUSDR&E 

Room  3D1079,  The  Pentagon 
Washington,  DC  20301 

Director 

Defense  Advanced  Research  Projects  Agency 
1400  Wilson  Boulevard 
Arlington,  VA  22209 


1 


1 

1 


1 


1 


1 


OUSDR&E 


1 


ATTN:  Dr.  G.  Gamota 

Deputy  Assistant  for  Research  (Research  in  Advanced  Technology) 
Room  3D1067,  The  Pentagon 
Washington,  DC  20301 

Director  of  Defense  Research  and  Engineering 
Engineering  Technology 
Washington,  DC  20301 


DIST-1 


mm 


DISTRIBUTION  (Continued) 


Director 

Defense  Advanced  Research  Projects  Agency/STO 
ATTN:  Commander  T.  F.  Wiener 
D.  W.  Walsh 
1400  Wilson  Boulevard 
Arlington,  VA  22209 

Commander 

US  Army  Aviation  Systems  Command 
12th  and  Spruce  Streets 
St.  Louis,  MO  63166 

Director 

US  Army  Air  Mobility  Research  &  Development  Laboratory 
Ames  Research  Center 
Moffett  Field,  CA  94035 

Commander 

US  Array  Electronics  Research  &  Development  Command 
ATTN:  AMSEL-TL-T,  Dr.  Jacobs 

AMLEW-R,  Henry  E.  Sonntag 
Fort  Monmouth,  NJ  07703 

Director 

US  Army  Night  Vision  Laboratory 
ATTN:  John  Johnson 
John  Deline 
Peter  VanAtta 
Fort  Belvoir,  VA  22060 

Commander 

US  Army  Picatinny  Arsenal 
Dover,  NJ  07801 

Commander 

US  Army  Harry  Diamond  Laboratories 
2800  Powder  Mill  Road 
Adelphi ,  MD  20783 

Commander 

US  Army  Foreign  Science  and  Technology  Center 

ATTN:  W.  S.  Alcott 

Federal  Office  Building 

220  7th  Street,  NE 

Charlottesville,  VA  22901 

Commander 

US  Army  Training  and  Doctrine  Command 
Fort  Monroe,  VA  22351 


No.  of 
Copies 


1 

1 


1 


1 


1 

1 


1 

1 

1 


1 


1 


1 


1 


DIST-2 


DISTRIBUTION  (Continued) 


No.  of 
Copies 


Director 

Ballistic  Missile  Defense  Advanced  Technology  Center 

ATTN:  ATC-D  1 

ATC-0  1 

ATC-R  1 

ATC-T  1 

P.  0.  Box  1500 
Huntsville,  A1  35808 

Commander  1 

US  Naval  Air  Systems  Command 
Missile  Guidance  and  Control  Branch 
Washington,  DC  20360 

Commander  1 

US  Naval  Air  Development  Center 
Warminster,  PA  18974 

Chief  of  Naval  Research  1 

Department  of  the  Navy 
Washington,  DC  20301 

Commander 

US  Naval  Ocean  Systems  Center 

Code  6003,  Dr.  Harper  Whitehouse  1 

San  Diego,  CA  92152 

Director 

Naval  Research  Laboratory 

ATTN:  Dave  Ringwolt  1 

Code  5570,  T.  Gialborinzl  1 

Washington,  DC  20390 

Commander 

Rome  Air  Development  Center 
US  Air  Force 

ATTN:  James  Wasielewski,  IRRC  1 

Griff iss  Air  Force  Base,  NY  13440 

Commander 
US  Air  Force 
AFOSR/NE 

ATTN:  Dr.  J.  A.  Neff  1 

Building  410,  Bolling  Air  Force  Base 
Washington,  DC  20332 


DIST-3 


DISTRIBUTION  (Continued) 


Commander 

US  Air  Force  Avionics  Laboratory 
ATTN:  D.  Rees 

W.  Schoonover 
Dr.  E.  Champaign 
Dr.  J.  Ryles 
Gale  Urban 
David  L.  Flannery 

Wright  Patterson  Air  Force  Base,  OH  45433 

Commander 

AFATL/LMT 

ATTN:  Charles  Warren 

Eglin  Air  Force  Base,  FL  32544 

Environmental  Research  Institute  of  Michigan 
Radar  and  Optics  Division 
ATTN:  Dr.  A.  Kozma 

Dr.  C.  C.  Aleksoff 
Juris  Upatnieks 
P.0.  Box  8618 
Ann  Arbor,  MI  41807 

IIT  Research  Institute 
ATTN:  GACIAC 
10  West  35th  Street 
Chicago,  IL  60616 

Dr.  J.  G.  Castle 
9801  San  Gabriel,  NE 
Albuquerque,  NM  87111 

Commande  r 

Center  for  Naval  Analyses 
ATTN:  Document  Control 
1401  Wilson  Boulevard 
Arlington,  Va  22209 

Dr.  J.  W.  Goodman 
Information  Systems  Laboratory 
Department  of  Electrical  Engineering 
Stanford  University 
Stanford,  CA  04305 

Eric  G.  Johnson,  Jr. 

National  Bureau  of  Standards 
325  S.  Broadway 
Boulder,  CO  80302 


DIST-4 


DISTRIBUTION  (Continued) 


Dr.  Nicholas  George 
The  Institute  of  Optics 
University  of  Rochester 
Rochester,  NY  14627 

Naval  Avionics  Facility 
Indianapolis,  IN  46218 

Dr.  David  Casasent 
Carnegie  Mellon  University 
Hamerschage  Hall,  Room  106 
Pittsburg,  PA  15213 

Professor  Anil  K.  Jain 
Department  of  Electrical  Engineering 
University  of  California,  Davis 
Davis,  CA  95616 

Terry  Turpin 

Department  of  Defense 

9800  Savage  Road 

Fort  George  G.  Meade,  MD  20775 

Dr.  Stuart  A.  Collins 
Electrical  Engineering  Department 
Ohio  State  Univeristy 
1320  Kennear  Road 
Columbus,  OH  43212 

US  Army  Materiel  System  Analysis  Activity 
ATTN:  AMXSY-MP 

Aberdeen  Proving  Ground,  MD  21005 

US  Army  Night  Vision  Laboratory 
ATTN:  DELNV-L ,  Dr.  R.  Buser 
Ft.  Belvoir ,  VA  22060 

Dr.  F.  T.  S.  Yu 

Penn  State  University 

Department  of  Electrical  Engineering 

University  Park,  PA  16802 

Dr.  William  P.  Bleha 

Liquid  Crystal  Light  Valve  Devices 

Hughes  Aircraft  Company 

6155  El  Camino 

Carlsbad,  CA  92008 


DIST-5 


V  -J-T-V  'f'K.7 


■"»  -  I  *VT m  "  '■  I  Jl  -W?VF!MP 


K 


V*. 

■w\ 

V*. 

1 


DISTRIBUTION  (Concluded) 


US  Army  Missile  Command 
AMCPM-PE-E,  Mr.  John  Pettitt 


-PE 

AMSMI-LP , 

Mr.  Voigt 

-0 

-Y 

-R, 

Dr.  M.  McCorkle 

-RN, 

Mr.  J.  Hagood 

-RE, 

Mr.  W.  Pittman 

-RD 

-RG 

-RG, 

Mr.  J.  A. 

McLean 

-RR, 

Dr.  R.  L. 

Hartman 

-RR, 

Dr •  J •  S • 

Bennett 

-RR, 

Dr  •  «J  •  G  • 

Duthie 

-RR, 

Dr.  D.  A. 

Gregory 

-RPR 

-RPT 

,  Record  Copy 

DIST-6 


No.  of 
Copies 


'VV-V1 


1 

1 

1 

1 

1 

1 

1 

1 

3 

1 

1 

1 

1 

1 

40 

15 

1 


.v 


•  *  * . 


-  V-* 

mULI 


■'.v 

WV.\ 
i>  .‘A 

*&*  ’’ 


5$ 

>  V 

\  V 

A  m\k 

-\A_V1 

—  » 


